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Abstract 

Perturbative QCD formula for inclusive energy spectra of heavy quarks from heavy quark 
initiated jets which takes into account collinear and/or soft logarithms in all orders, the 
exact first order result and two-loop effects is applied to distributions of heavy flavoured 
hadrons in the framework of the LPHD concept. 
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1 Introduction 



Heavy flavour physics is now extensively studied experimentally at both e + e~ and hadronic col- 
liders. Experiments at Z° have led to the availability of new data on the profiles of jets initiated 
by heavy quarks lira. Further progress is expected from the measurements at LEP 2 and, espe- 
cially, at a future linear e + e~ collider. The principal physics issues of these studies are related not 
only to testing the fundamental aspects of QCD, but also to their large potential importance for 
measurements of heavy particle properties: lifetimes, spatial oscillations of flavour, searching for 
CP-violating effects in their decays etc. Properties of b-initiated jets are of primary importance 
for analysis of the final state structure in ti production processes. A detailed knowledge of the 
b-jet profile is also essential for the Higgs search strategy. 

The physics of heavy quarks has always been considered as one of the best testing grounds 
for QCD. Despite this, not many attempts to derive self-consistent perturbative (PT) results 
for the profiles of heavy quark jets have appeared so far. This paper is aimed at a condensed 
presentation of the part of the results of the long-term Leningrad/St. Petersburg QCD group 
project^ that concerns PT analysis of the inclusive energy spectra of heavy quarks. 

The selection of this subject is motivated by the following topical questions: 

i) To what extent can heavy quark distribution be treated as a purely perturbative (Infrared 
Safe) quantity? 

ii) Can we describe the gross features of jets initiated by heavy quarks without invoking 
phenomenological fragmentation schemes? 

iii) What kind of information can be extracted from the measurements of heavy quark energy 
spectra and, in particular, from the scaling violation effects? 

There are two ingredients of the standard Renormalization Group (RG) approach to the descrip- 
tion of the energy spectra of heavy flavoured hadrons Hq. Here one starts from a phenomeno- 
logical fragmentation function for the 

Q(x) -> H Q {x H ) (1.1) 

transition and then traces its evolution with the annihilation energy W by means of PT QCD. 
Realistic fragmentation functions H exhibit a parton-model-motivated maximum at (hereafter 
M is the heavy quark mass) 

xh const 
1 ~ . 

x M 

Taking gluon radiation at the PT-stage of evolution into account would then induce scaling 
violations that soften the hadron spectrum by broadening (and damping) the original maximum 
and shifting its position to larger values of 1 — xh with W increasing. 

This approach has been successfully tried by B. Mele and P. Nason @, who have studied 
the effects of multiple soft gluon radiation and have been looking for realistic fragmentation 
functions to describe the present day situation with heavy particle spectra, and to make reliable 

1 for some basic references see || || 
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predictions for the future. Being formally well justified, such an approach, however, basically 
disregards the effects that the finite quark mass produces on the accompanying QCD brems- 
strahlung pattern, since the W-evolution by itself is insensitive to M (save the power suppressed 
M 2 /W 2 corrections). 

At the same time, as far as one may consider M a sufficiently large momentum scale, it is 
tempting to carry out a program of deriving the predictions that would keep under PT control, 
as much as possible, the dependence of Hq distributions on the quark mass. 

In Section 2 we present and discuss the PT formula for inclusive quark spectra, D(x; W, M), 
that in the QED case would give an unambiguous all-order-improved absolute prediction for the 
inclusive energy distribution of muons produced in 



Given an explicit dependence on m M , one may theoretically compare, say, the /i and r spectra, a 
goal that, generally speaking, can not be achieved within the standard "evolutionary" framework. 

In Section 3 the topical questions listed above are addressed. We discuss the problem of 
"infrared instability" of the PT quark spectra and suggest a solution of this formal difficulty 
based upon the notion of an "infrared regular" effective coupling. Having regularized the PT 
expression in this way we observe that the gluon bremsstrahlung effects (Sudakov suppression 
of the quasi-elastic kinematics x — *■ 1) lead to particle spectra that have a similar shape to the 
non-PT fragmentation function. 

An emphasis of the role of PT dynamics has been successfully tested in studies of light 
hadron distributions in QCD jets. Prompted by the success of the MLLA-LPHD approach EJ, 
we consider an option of describing the gross features of heavy quark initiated jets entirely by 
means of PT QCD, without invoking any fragmentation hypothesis. 

The LPHD philosophy would encourage us to continue the PT description of an inclusive 
quantity down to small momentum scales, and with the hope that such an extrapolation of the 
quark-gluon language would be dual to the sum over all possible hadronic excitations. For the 
problem under consideration, such an approach could describe the energy fraction distribution 
averaged over heavy-flavoured hadron states, the mixture that naturally appears, e.g., in the 
study of inclusive hard leptons. 

In Section 4, we report results on the comparison of the PT predictions with the measured 
c and b energy losses at different annihilation energies. Fitting the rate of scaling violation in 
(x) (W) results in the QCD scale parameter A that, after being translated into the popular MS 
scheme, agrees with the values extracted from other studies. On the other hand, the absolute 
values of (x Cj b) at a given energy are most sensitive to the behaviour of the radiation intensity, 
say, below l-2GeV. 

Limited experimental information does not at present allow us to disentangle various possible 
shapes of af near the origin. At the same time, a quite substantial difference between (x c ) and 
(xb) fits nicely into PT-controlled mass dependence, provided the characteristic integral of the 
effective coupling over the "confinement" region assumes a fixed value, 



e + e — > ^(x) /i T + . . . 



(1.2) 
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(with k the linear momentum variable). Given this value, one may predict the differential c and 
b energy distributions at arbitrary W, predictions which prove to be unaffected in practice by 
our ignorance of the detailed behaviour of a e f in the small momentum region. 

In Section 5 the results of the study are discussed and conclusions drawn. 

Detailed derivation of the master formula for the inclusive heavy quark energy spectrum is 
given in Appendices. In Appendix A the problem of running coupling in the perturbative radiator 
is dealt with. Appendix B is devoted to the second loop effects in the anomalous dimension. 



2 Perturbative Energy Spectrum of Leading Quarks 

We denote the cms annihilation energy by W, and M = m ■ W is the quark mass. 
Here we describe perturbative QCD result for the inclusive energy spectrum of heavy quarks, 
D(x; W, M), x = 2Eq/W, produced in e + e~ — > Q(x) + Q+ light partons. 
This function has the following properties: 

• it embodies the exact first order result = a s ■ f(x; m) and, as a consequence, 

• it has the correct threshold behaviour at 1 — 2m <C 1; 

• in the relativistic limit m <C 1, it accounts for all significant logarithmically enhanced 
contributions in high orders, including 

1. running coupling effects, 

2. the two-loop anomalous dimension and 

3. the proper coefficient function with exponentiated Sudakov-type logs, which are es- 
sential in the quasi-elastic kinematics, (1— x) <C 1; 

• it takes into full account the controllable dependence on the heavy quark mass that makes 
possible a comparison between the spectra of b and (directly produced) c quarks. 

By "threshold behaviour" we mean here the kinematical region of non-relativistic quarks, \W — 2M\ 
M, in which gluon bremsstrahlung acquires additional dipole suppression. At the same time we 
will not account for the Coulomb effects that would essentially modify the production cross 
section near the actual threshold, W — 2M <C M 

We are mainly interested in the "leading" energy region of x of the order of (or close to) 
unity. For this reason we disregard in what follows potentially copious production of QQ pairs 
via secondary gluon splitting^, i.e. singlet sea contribution, as well as the specific Q — > Q + QQ 
transition that appears in the non-singlet anomalous dimension beyond the first loop (for a 
review see 



15| and references therein). 



Therefore the quark distribution we define below will satisfy the sum rule 

1 dxD(x;W,M) = 1 . (2.1) 

2m 



2 apart from the integral effect of the unregistered pairs embodied into the running coupling, see below 
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The accuracy of the PT prediction described in this paper is 

D(x;W,M) ■ [l + 0(a£) 



(2.2) 



This expansion is uniform in x, that is, the a s correction does not blow up (neither as a power, 
nor logarithmically) when (1—x) — ► 0. 



2.1 Definition of the leading quark distribution. 

The inclusive cross section for single particle production in e + e~ annihilation is conveniently 
written in terms of structure functions as, e.g., @ 



_i da 



(2.3) 



with a the standard Born cross section factor. 

We define the normalized inclusive energy spectrum as 

_ 1 da _ f dj 



(2.4) 



where the contour T runs parallel to the imaginary axis in the complex moment j plane. 

Aiming at an exact account of the first order QCD effects, we include in crtot the first a s 
correction to the annihilation cross section, that in the case of relativistic quarks reads 



0tot = °o 



l+*C F a + Ofa 2 )] ; a = a(W 2 ) = 



a s (W) 



-i 



2tt 2N c 3 

With this definition, the structure functions JF in the moment representation become 

^(j) = c Q (j) (l-fCVafr 1 -1]) D r [l + (a 2 ) ] , 
clU) = Cpaj" 1 , c T (j) = 1 + \Cpaj~ 1 , c 2 (j) = 1 - \C F aj 

2.2 The radiator. 

The integrand of ( |2.4j ) may be written in the exponential form as 

dw(x; W, M) 



(2.5) 



lnZ^ 





' x j - 1 - 




/ dx 


- 1 


hm 







dx 



(2.6) 



where the "radiator" dw/dx originates from the improved first order gluon emission probability, in 
which finite mass and running coupling effects have been included^. The corresponding expression 
reads 



C~ l v — 
F dx 



Q 2 dt 

K 2 t 



a(t) 
2x 



2(x-2m 2 ) 
1 — x 



+ (-\l-x) 



a'(t)(l-x) + a 2 (t)A^ 2 \x) (2.7a) 



1—x 



a(Q 2 ) + a(n 2 ) +C 



„! x (x — 2m 2 ) 



2(1 — x) VI — x+m 2 



1 — x 



a(Q 2 



(2.7b) 



'Detailed derivation of the radiator is presented in Appendix A; see also 
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Here the two characteristic momentum scales have been introduced, 

Q 2 =Q 2 (x) = W 2 - {1 ~f 2 zo, 
1 — x + m z 

k 2 =k 2 {x) =M 2 - {l ~ X)2 , 

with 

zq = | (x — 2m 2 + \J x 2 — Am 2 J = a; + O (m 2 ^ 
The following notations were also used: 
M 



(2.8) 



m 



W 



< \ ; u = Vl - 4m 2 ; /3(a;) = a/1 - Am 2 /x 2 , < /? < v ; C = 1 + 2m 2 



cilnfc 2 v " v y 2vr 

The radiator (12.71) vanishes at 



(2.9) 

x = x min = 2m; ( /3(z min ) = , Q 2 = H/M (1 ~ 2m) = k 2 ) , 
y 1 — m y 

thus justifying the lower kinematical limit in eq. ( |2.6|) . 

Relativistic Approximation. In the relativistic approximation, m <C 1, we set v — (3 — Q — 1 
to get a simplified expression 

C ^ = if f { fl(t)P(x) " a ' W (1_a:) + fl2 W A(2) ^) } 

f-2r r 2 1 r-2r^ (2.10a) 

+ a(Q 2 ) -^ + -^_ + a(K* ' " 



Q 2 = W 2 x(l - x) , k 2 = M 2 (l — x) 2 / x , (2.10b) 

where 

1 4- r 2 

P(z) = -L-. (2.11) 

The integration variable t determines the physical hardness scale of the running coupling, and is 
related to the transverse momentum of the radiationF]. In the dominant integration region, 

k 2 ± ^W 2 , t = x-k\. (2.12) 

The lower limit t > k 2 sets the boundary for the essential gluon emission angles, 

t ~ E Q ^ ■ (^6) 2 > k 2 ~ M 2 |^ -| =► 6 > M/E = O . 

This restriction manifests the "dead cone" phenomenon characteristic for bremsstrahlung off a 
massive particle. It is largely responsible for the differences between radiative particle production 
in jets produced by a light and a heavy quark (excluding the decay products of the latter), see 

!_!■ 

4 For the space-like evolution similar relation holds with a; -1 substituted for x in ( |2. 12 ) 
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2.3 Logs and their exponentiation. 

Singularities of the radiator ( |2.7|) at x — 1, when driven through the inverse Mellin transform 
Q2.4|) , ( |2.6|) , give rise to ln(l— x) terms. Bearing this in mind, one may represent the term-by-term 
structure of the radiator by the following symbolic expression, 

(2.7a) =^ a{CS + S 2 ) + aC + a 2 C + [ a 2 C + a 2 {CS + S 2 ) ] ; (2.13a) 
aS + aS + aS , (2.13b) 

where C = \nW/M and S = In 1/(1— a;) represent large logs that usually reflect enhancements 
due to quasi-collinear and soft gluon radiation respectively. As we shall see shortly, under the 
"physical" definition of the QCD coupling the last 2nd order term in ( |2.13a| ) is free from the 



double-log contribution (the underlined piece) so that the "convergence" of the PT expansion 
( [2.13j ) is improved. 

Exponentiation of collinear logs follows from the general factorization theoremF]. Arguments 
in favour of exponentiation of aS 2 + aS contributions according to ( |2.6| ) have been given in 
|/J, motivated by the factorization property of soft bremsstrahlung . Exponentiation of non- 
logarithmic O (a) corrections, corresponding to the terms in the non-integral part of the radiator 
(|2.7b ) that are regular at (1— x) — > 0, is questionable (indeed wrong) and unnecessary within the 



accuracy adopted (see (|2.2D). Keeping such terms in the exponent is a matter of choice. We do 
so just to simplify the result and make the normalization sum rule (|2.1|) automatically satisfied, 



When translated into the standard language of the RG motivated approach, the integral part 
of the radiator ( p.lOj ) may be said to embody the anomalous dimension together with (a part of) 
the correction to the hard cross section due to the x-dependent factor in the upper integration 
limit: 

r W 2 x{i-x) A+ f w 2 A+ /r7 2 ln 2 n — t)\ 
J ja{t)P{x)^J ja{t)P{x) + a{W 2 )P{x)ln[x{l-x)]+0 K _ M . (2.14) 

Non-uniformity of such an expansion with respect to the potentially large logarithm, ln(l — x), 
explains why we preserve the exact kinematical limits in ( |2.7| ). By neglecting the last term in 
( p,14j ), one would lose control of the a 2 S 3 terms. Such contributions formally belong to the 2nd 
order correction to the hard cross section, and therefore lie beyond the reach of the two-loop RG 
analysis. In the meantime, such an ignorance would undermine the possibility of keeping track 
of essential first subleading corrections, a n log 2n_1 , at the level of running coupling effects in the 
quark form factor (Sudakov suppression). 

Let us stress that eqs.( |2.6p , ( |2.7| ) solve the problem of all-order resummation of soft radiation 
effects both in the hard cross section^ and in the coefficient function, that is at low momentum 
scales K 2 oc M 2 . In particular, the term in ( 2.7b|) proportional to a(n 2 ) (which is ^-independent 



and therefore gets lost in the standard evolution approach) accounts for the exponentiated contri- 
bution from soft gluons at small emission angles < 0q, k\ M 2 — the dead cone subtraction 



5 this applies to C-contaminated terms of ( 2.13a ) as well as to the last contribution in ( 2.13b ), that actually is 
due to hard gluons collinear to the Q momentum (the "backward jet" correction 01) 
6 for a review of such resummation programs see, e.g., [T^ | and references therein 
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effect. Notice that the similar term proportional to a{Q 2 ) (which does belong to the hard cross 
section corrections) is due to the dead cone subtraction in the backward jet. 



2.4 Second loop effects. 

The preliminary result of []7J is improved by taking into full account effects of the two-loop 
anomalous dimension embodied in the a' and A^ 2 ^ of ( p.7a|) . We prefer to treat these two terms 
separately, as the former naturally emerges in the dispersion relation approach to definition of the 
running coupling @ . It shows that beyond the first loop the "soft" and the "hard" parts of the 
GLAP splitting function ( 2.11 ) in the anomalous dimension actually acquire different physical 
interaction strength @, 
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(i) 



1 + x 2 
1 — X 
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(2) 



2x 
1 — x 



+ [a 



x) 



The actual expression for A depends on the renormalization scheme chosen for the one-loop 
coupling. It may be fixed, e.g., by matching with the known two- loop MS result for the non- 
singlet fragmentation function. To this end one has to consider the VT-dependent part of the 
relativistic radiator (|2.10|) , 



C 



+ a(W 2 ) (f^ + ) + const + O (a 2 (W 2 

{ 1 — x 2(1— x) J v 



(2.15) 



and, with account of (|2.5|) relating D to the structure functions, compare the scaling violation 
rate according to Q2.6|), fl2.15| ) with that computed in [TP|-f2"T|]. 

By doing so (see Appendix B) one arrives at (Ca = N c = 3 , T R = |) 



AE(x1 



P{x)-K + A (2) ( 



x) 



A (2) i 



K 



X) 



67 
18 



7T" 

~6 



10 



n f Tj 



R 



( 4.565 
4.010 
3.454 



(n f = 3) , 
(n f = 4) , 
{n f = 5) ; 



C F V{x) + K(x) 



(2.16a) 

(2.16b) 
(2.16c) 



V term is responsible for the difference between the time-like and space-like anomalous dimen- 
sions. Explicit expressions for V(x) and TZ(x) are given in Appendix B. 

(2) 



Introducing the Physical Coupling. In the x —> 1 limit, A^ (|2.16| ) peaks together with 



P(x) oc (1 — x) 1 . At the same time, A'- 2 -' is less singular in this limit, 

V(x)/A (2) (x) oc (1 -x) ln(l 
TZ(x)/A (2 \x) oc (1-x) 2 . 



x 



7 



As far as large quark energies are concerned, In 1/(1— x) > 1, the A^ 2 ^ term constitutes a small 
correction. Thus, the major part of the two-loop effect according to ( j2.16Q reduces to a finite 
renormalization of the leading term aP(x). This naturally suggests absorbing the correction by 
introducing the "physical" effective coupling according to 



P(x) + a 2 A^(x) = a cft P(x) + a 2 A (2 \x) + £> (a 3 ) ; (2.17a) 
a Ms( 1 + a ^)' a Ms ~ a <=ff (1 ~ a cff ■ (2.17b) 



Eq.(^TTTJ) relates the scale parameters A of the two schemes. This relation within two- loop 



accuracy reads 




lnAws + — = hi Arm + 



MS ' a ~ 1±1J1 MS ~ 11 r< 4 ^ r j i 



^ u f 1 R 

3 } (2.18) 

1.66 (n/ = 3), 
1.57 (n/ = 5) . 

Our A eff coincides with the A^q introduced by Catani, Marchesini and Webber in |2^ . 



Quark Thresholds in the Running Coupling. One additional comment is in order con- 
cerning the essence of a c f . In spite of the formal equivalence of the two representations (|2.17a|) , 



the former one is physically preferable since A^(x) is free from an ill-defined quantity rif repre- 
senting the number of "active" quark flavours. The t-integration in (|2.7|) runs over a broad region 
that is sliced by the finite quark mass scales. Naturally, one has to increment rif when passing 
such a scale, e.g. nf — 3 — > rif — 4 when going through t ~ 1.5GeV. A^jg has to be redefined by 
adjusting to a new /5-function value, and Aj^g should be changed accordingly. The problem ap- 
pears to be entirely technical, however, and reflects a deficiency of the MS prescription as one of 
the schemes based upon the dimensional regularization technique that eventually treats fermions 
as massless particles: The discontinuous behaviour of A and A gets compensated in the physical 
combination (|2 . 1 7a|) that emerges in the full two-loop anomalous dimension. It is worthwhile to 
notice that our convention (rif independent A^(x), the scale parameter A defined by Q2.18 )) is 
in accord with the BLM prescription @ for optimizing the choice of "physical coupling" . 

For the argument sake, one may once again invoke the QED example ( |1.2|) . In this context, an 
application of the MS scheme would reveal the same problem with the lepton thresholds (e, /z, r 
etc. ), while a cB would be nothing but the physical running QED coupling, the one given by the 
Euclidean photon renormalization function Z$(t) that is unambiguously linked by the dispersion 
relation to fermion pair production (with the finite fermion mass effects fully included). 

The improved one-loop expression for a eS accounting for finite quark mass effects reads (for 
a detailed discussion see HT6f ) 



o.r'(Q 2 ) = ft In ^ -^Enll + ^M + a .«rV), 0.«A/,); (2.19a) 

n « WM?) r l. (2.19b) 
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where j3e accounts for gluon and massless quark contributions (j3g = &[ n/ =s] = 9) and II is the 
standard fermion loop polarization operator known from the mid-fifties (see below eq. ( |3.24j )). 
According to ( 2.19b|) , heavy flavours decouple at low momentum scales (as they ought to). In 



the "ultraviolet" regime, II acquires a constant subtraction from the leading log behaviour, which 
term, if treated as the two- loop correction, explains the rif dependent piece of the /C factor ( |2.16b|) 
that enters the relation ( |2.17b| ) between the physical and the MS couplings: —^T R = |Tr- (— |). 

A^ 2 ) is Practically Negligible. As an estimate of the magnitude of the correction we 
present its contribution to the second moment j = 2 that describes quark energy losses: 

C F {a m Pj(l + aJC) + a 2 hf) = C F Pj{a cS + a 2 6j} + O (a 3 ) ; 

31 (2-20) 
C F + — (C A -2C F ) = -0.1735. 
3o 



21 /85 _ vr 2 * 
36 + \ 36 ~ y 



Being numerically very small already for j — 2, the relative correction 5j continues to fall as j 1 
with j increasing. Therefore in practice one may use (|2.7|) with A^ 2 ) = 0, provided the physical 
effective coupling ( |2.17p| ) is used. 



2.5 Axial vs. Vector currents in QQ production. 



The master equation ( |2.7| ) has been written for the case of QQ production via the vector current, 
in which case 

( = (vS fEa^ = 1 + 2m . = . (3 _^. 

This factor tends to 1 in the relativistic limit m — > but depends otherwise on the production 
channel (see Appendix A for details). It is worthwhile to note that beyond the leading twist 
approximation, generally speaking, there is no way to define the universal fragmentation distri- 
butions that would be independent of the production stage. The origin of such non-universality 
(that shows up at the level of a s m 2 In m 2 terms) may be traced back to the process dependent 
radiation of hard gluons,[] 

da oc ujg ■ dujg oc (1 — x) dx. 

In practice we have to deal with a mixture of vector and axial production currents. In the pure 
axial case ( should be substituted by 

a = ^T = l-4m 2 = , 2 . 

ff [A](l) 

In addition, the C _1 (l ~ x) term in (|2.7|) acquires an extra mass correction factor (l+2m 2 ). Both 
this effect and the difference between (y and (a are proportional to m 2 = (M/W) 2 . We conclude 
that 



7 both the "soft" and "semi-soft" terms of the radiation spectrum, oc ui g 1 • dio g and oc 1 • du> g , remain 
universal 0@. 
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• the difference between the radiation spectra in V- and A-channels vanishes as (1 — v 2 ) in 
the relativistic case, while 

• for non-relativistic quarks, »<1, the axial contribution to the cross section is relatively 
suppressed as A/V ~ v 2 . 

Therefore (|2.4| )— fl2~77|) may be used in practice for the realistic V + A mixture, in particular, at 
the Z° peak. 



3 Hadronization and Infrared Finite a, 



In the QCD context the non-PT phenomena inevitably enter the game. It is important to stress 
however that the strong interaction not only determines the hadronization transition ( |1 . 1| ) but 
also affects to some extent the quark evolution stage: The high order effects embodied into the 
running coupling seem to undermine the very possibility of the PT analysis. 

It is often believed that the large quark mass, M>A, provides a natural cut-off, which keeps 
the relevant space-time region compact enough to avoid the truly strong, non-PT interaction in 
a course of the quark evolution. 

To illustrate the point we invoke a rough estimate that follows immediately from (|2.6|) and 
is valid for numerically large x values in the double-log approximation , 



D(x;W,M) oc (1-x) 



with A£ the characteristic evolution integral 



(l-x) 2 M 2 kj_ 



71 



In the fixed coupling approximation we would get 



D(x) oc (1 



x) 



■ exp 



C F a s 



2tt 



\n 2 (l-x 



(3.1a) 



(3.1b) 



(3.2) 



In the small coupling regime, Cf— In < 1, this distribution exhibits a sharp peak at large x 
followed by a steep fall off in the x — > 1 limit due to the double logarithmic Sudakov suppression. 
Qualitatively similar shape for the heavy hadron spectrum due to Q — > Hq transition one expects 
from the parton model considerations cJ which have lead to the so called Peterson fragmentation 
function o 



Cq(x) 



N 

x 



1 — x 



-2 



x 1 — x 

This function peaks at 1 — x ~ e 1//2 <^ 1 leading to 



N 



4e 



1/2 



7T 



l + O 6 



.1/2 



(1-X) 



,1/2 



fragm 



(lne^-l) [l+0(e 



,1/2 



oc e 



1/2 



(3.3) 



(3.4) 
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3.1 Peterson fragmentation function vs. integrated coupling. 



One may single out effects of the non-PT momentum region by factorizing the D spectrum in 
the moment representation ( |2.6|) into the product of the "safe" PT part and the "confinement" 
factor, 



Da = Di 



x dP 



(3.5a) 



In other words, we split the radiator into two pieces corresponding to large and small transverse 
momentum regions (t w k\ for x close to 1), 



dw 
dx 



dw 
dx 



k 2 ± >fi 2 



+ 



dw 



dx } 



(C) 



k\<^ 



(3.5b) 



This formal separation becomes informative if one is allowed to choose the boundary value \i well 
below the quark mass scale (e.g., \i— lGeV providing /i/M<^C 1 for the b quark case). Within 
such a choice only x close to 1 would contribute to w^ c \ For illustrative purposes let us neglect 
subleading fi/M effects and retain the most singular term only to get an estimate 



dw) iC) 
dx] 



JL 

M 



[l-x. 



dt , , 2C p 
a (t) —L. + 

[(l-x)M] 2 t 1 — X 



(3.6a) 



One arrives at 



ln-D 



(C) 



dz 



y- 1 



1 



f 1 dt , , 2C t 
— a(t) — - 

zM] 2 t Z 



2C f 



t* dk a s (k) f k du 
o k 



71 



u 



—) 



3-1 



(3.6b) 



It is important to stress that it is the first "confinement insensitive" factor of (|3.5a|) only that 
depends on W. Therefore, the ratio of the moments 



Dj(W, M)/ Dj(W Q , M) 



(3.7a) 



as a function of W and the heavy quark mass is expected to be an "infrared stable" PT prediction. 
Another message one receives observing the structure of the radiator (|2.7|) is that the ratio of 
the moments for different quarks, 



D j (W,M 1 )/D j (W,M 2 ) 



(3.7b) 



should tend to a ^-independent confinement sensitive (sic\) constant in the relativistic limit 

We proceed with the estimate of the "confinement" factor ( |3.6b| ). As far as ///M may be 
treated as a small parameter, for finite moments j ~ 1, jfi/M 1, 



v dk a s (k) 



71 



1 + 



-2C F (j - 1) ^ , (3.f 
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with 

r*» a.(k) 



H a = / dfc . (3.9) 

JO 7T 



In particular, for the energy losses, j = 2, one has 

1 + 



lni}f>=ln<x>^=-2<^ 
2 \ / M 



M 



(x) (C) =exp{-2C F ^j. (3.10) 



If ( p.8|) were applicable for all j, the inverse Mellin transform would result in a singular distri- 
bution 

D{xi c) = 5 (x - (xf) . (3.11) 

This singularity gets smeared when a proper treatment is given to the large j region. To this 
end one may use a simple expression that interpolates between ( |3.8|) and the correct logarithmic 
asymptote of the "confinement radiator" in ( |3.6| ), 



InDf «-2<^r^^ffiln 
J Jo k n 



k t \ 



(3.12) 



As a result, a distribution emerges that peaks around x = (x}^ and is rather similar in shape to 
( [3. 3D - The extreme x— ► 1 asymptote is determined by the region of parametrically large moments 
(j) oc (1 — x)^ 1 ^> M/fi. It is different for two regimes: 



a s (0)>0, InDf ~ -O^i In" =► £(x) oc ^ ^;(3.13a) 

J 7T V MJ 1 — X 

« s (0) = 0, lnD^~-2C^Eoln(^) D(z) oc (1 - x) 2 ^ 3 "- 1 , (3.13b) 
where, in the latter case, 



^ dfc a s (k) .„ _ 

^ = i < oo. (3.14) 

JO K 7T 

Reproducing the concrete behaviour of the Peterson function ( |3.3|) in the large x limit, C(x) ~ 
(1 — x) 2 , would require S = | in the second regime ( |3.13b| ). We conclude that the particle 
distribution originating from (|3.6|) is capable of reproducing the gross features of the popular 
Peterson fragmentation function (provided, naturally, that the notion of the infrared regular 
effective coupling is implanted in the PT- motivated "confinement" radiator). 

An important message comes from comparing the energy losses that occur at the hadroniza- 
tion stage. Based on the pick-up hadronization picture, the characteristic parameter e in ( |3.3| ) 
has been predicted to scale @ as 

6Q * (m) 2(xM " 2 ' (3 - 15) 

(with m q the quantity of the order of constituent light quark mass.) Confronting ( |3.4j ) that (up 
to a logarithmic factor) scales as - /6q with the PT prediction ( |3.10D we get 



(1 " x) fragm - <=► l-(x) {C) = l- exp {~2C F ^} * 2C F ^ 
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which justifies the expected scaling lawf[ 

It is worthwhile to remember that neither the Peterson function nor our PT-motivated "con- 
finement" distribution is an unambiguously defined object. The former as an "input" for the 
evolution is by itself contaminated by gluon radiation effects at the hard scale t ~ M 2 that 
are present even at moderate W > 2M. On the other hand, D\k\ < fi 2 } crucially depends 
on an arbitrarily introduced separation scale \i that disappears only in the product of the fac- 
tors responsible for "PT" and "non-PT" stages fl3.5a| ). Nevertheless, bearing this in mind, one 
may still speak of a direct correspondence between these two quantities, namely, C{x) in the j 
representation and as given by ( |3.6| ). 

This means that instead of convoluting phenomenological C (x) with the jy-dependent "safe" 
evolutionary quark distribution one may try to use consistently the pure PT description that 
would place no artificial separator between the two stages of the hadroproduction. From the 
first sight, one gains not much profit substituting one non-PT object — the phenomenological 
fragmentation function C(x) - - by another unknown, namely, the behaviour of the effective 
long-distance interaction strength a s (k) (at, say, k < 2GeV). There is, however, an important 
physical difference between the two approaches: a s should be looked upon as an universal process 
independent quantity. Therefore quite substantial differences between inclusive spectra of c- and 
b- flavoured hadrons should be under complete control according to the explicit quark mass 
dependence embodied in the PT formulae. 



3.2 Modeling the coupling. 

To study infrared sensitivity of PT results one can try different shapes of the effective coupling 
or, equivalently, different ways to extrapolate the characteristic function 

im = f £ «&> + const (3.16) 
to the "confinement" region of small Q 2 . Using the one-loop expression for the coupling, 

for £ one gets 

4 k 2 

f « (k 2 ) = - In In — + const , (3.18) 
b A 2 

which expression is defined only for k > A. 

For the two-loop effective coupling we use the standard approximate expression 

af\k) = of>(*) (l - 6lh ^f A2) of (*)) (3.19a) 



8 Similar behaviour was advocated recently by R.L. Jaffe and L. Randall who have exploited the difference 
between the hadron and the heavy quark masses as a small expansion parameter. 
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with 



h = y N 2 - (H]V c + 2C F ) n 



f 



(3.19b) 



and given by the one- loop formula (|3.17|) . Analytic expression for £ then reads 

k 2 



e (2) (k 2 ) 



- lnL + -— —1 + const 



L = In 



A 2 



(3.20) 



3.2.1 F-model. 

The simplest prescription which we refer below as the F-model consists of freezing the running 
coupling near the origin. One follows the basic PT dependence given by either (|3.17|) or (|3.19a|) 
down to a certain point k 2 where the coupling reaches a given value 



a s (k c 



A 



IT 



(3.21a) 



and then keeps this value down to k 2 = 0. £ then takes the form 

= i{k 2 ) ~ m) , k 2 > k 2 c ■ 
i(k 2 )=A\n(k 2 /k 2 c ) , k 2 <k t 

with k c related to A by ( |3.21a| ). 



2 

c ' 



(3.21b) 



3.2.2 G-model. 

The set of G p -models (Generalized shift models) gives another more flexible example for the trial 
effective coupling. It emerges when one regularizes the evolution function ( |3.18| ) as follows, 



e ] {k 2 



-In In [-— + C P 

b \A 2 p p 



+ const 



C P >1 



(3.22a) 



which corresponds to the effective coupling 



aW(Jfe) _ d^{k 2 ) 



TV 



dink 2 



k 2 P 



k 2 P + C v A 2 P 



P 



b In (k 2 P/A 2 P + C p 



(3.22b) 



This expression preserves the perturbative asymptotic form ( |3.17 ) up to power corrections 
A 2p /Q 2p . Notice that the effective coupling ( |3.22b[ ) with C p = 1 has a finite limit a s (0)/ir = Ap/b, 
while for C p > 1 it vanishes in the origin. 

For the two-loop coupling one substitutes in ( |3.2U| ) 



v 1 



(3.23a) 
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which results in 



a&{k) _ d^(k 2 ) 



71 



din k 2 



k 2p 



k 2 P + C v A 2 P 



4p ( bilnLp 



bL 



b 2 L r 



(3.23b) 



Trial shapes of the effective coupling in the G2 model (G-model with p = 2 and the two-loop a s 
with rif — 5 massless flavours) are displayed in Fig.|I] for different values of the parameter C<i- 
Crosses mark the curve that provides the best fit to mean energy losses (see below). 



G 3 — model Effective Coupling 



2.0 



n 1 1 1 1 1 1 r 



n 1 1 1 1 1 1 r 



a s eff (k) 



1.5 



1.0 



0.5 



0.0 



C 2 = 2 



A (5) = 580 MeV 




2.5 



Figure 1: Trial shapes of af in the G2 model (two-loop, rif = 5). Marked by crosses is the best-fit 
coupling. For comparison the best-fit F-model coupling (^4=0.19) is shown with dashed line. 



3.2.3 Quark Thresholds. 

When it comes to an accurate account of heavy quark thresholds in a s we modify the logarithmic 
denominator of a^\k 2 ) in ( |3.f 9a|) according to ( |2.19| ) as follows 

k 2 



6 M ln A2 ~ 3 



n 1 + — -± + n 1 + 



k 2 



k 2 



with 



„, 9 . v(3 — v 2 ) , 1 + v 9 
U(v 2 ) = ^ - ln - — - + v 2 



(3.24a) 



(3.24b) 
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Given infrared regular behaviour of a e f, numerical evaluation of the inverse Mellin transform 
with the PT radiator ( |2.7| ) becomes straightforward. Inclusive energy heavy quark spectrum 
obtained along these lines is concentrated (has a sharp maximum) near xq = 1 at W > 2M and 
softens due to gluon bremsstrahlung effects with W increasing. Similar softening one achieves 
increasing radiation intensity in the PT domain (by taking a larger value of A) and/or at small 
momentum scales (by varying "confinement parameter" of the model; larger A, smaller C p ). 

Since available experimental information on differential heavy quark spectra is rather scarce 
at present (and possibly contradictory), we restrict ourselves by considering the mean energy 
losses which, as we have discussed above, can be studied to quantify the influence of non-PT 
effects. 

The value of (xq) = Dz(W, M) shows up quite a strong variation with A/C. However, from 
the general factorization argument (|3.7a| ) one would expect the energy dependence of (xq) to 



stay well under PT control. As demonstrated in [|3^], the normalized quantity 

(x Q )(W)/(x Q )(W ) (3.25) 

is indeed practically insensitive to the variation of the "confinement parameter" of the model. 
Therefore the ratio ( |3.25| ) may be looked upon as an infrared stable quantity suitable for mea- 
suring A as has been suggested by P.Mattig 0. 

Position of the peak in the energy distribution seems to give another less trivial example 
of a stable prediction. Once again, as in the case of (xq), the absolute value of the peak po- 
sition depends strongly on the chosen A value. At the same time the normalized quantity 
x peak (W)/ x peak (Wo) exhibits much weaker A/C-dependence than the PT-controlled dependence 
on A f§. 

Thus the H^-evolution (scaling violation) in quark energy losses (xq) (W) allows one to extract 
the scale parameter A which value proves to be practically insensitive to the adopted scheme of 
af extrapolation. At the same time the absolute values of (xq) are quite sensitive to the gross 
effective radiation intensity below l-2GeV which makes it possible to quantify corresponding 
"confinement parameter" of the scheme and, thus, the shape of the coupling. 

Worthwhile to notice that there is a natural theoretical scale the "measurement" of af below 
l-2GeV to be compared to. As shown by Gribov HI, in the presence of light quarks colour 
confinement occurs when the effective coupling (parameter A of the F-model) exceeds rather 
small critical value 

f n ^ crit r -i 

A > pj =Cp l 1-yf w0.14. (3.26) 

Thus, within the Gribov confinement scenario an interesting possibility arises. Namely, if phe- 
nomenological af extracted from the data does exceed af lt but remains numerically small this 
would provide a better understanding of the PT approach to multiple hadropro duct ion in hard 
processes. 



4 Numerical Analysis of Energy Losses 

In this Section we compare experimental data with the generalized PT prediction which embodies 
the notion of the infrared regular effective QCD coupling. As an input we take the world average 
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values of (xq) listed in Table [p. 



process 


^ cm (GeV) 




1. 






10.4 


0.727 ±0.014 


2. 


c - 


-+ D* + ... 


30 


0.587 ±0.015 


3. 






91 


0.508 ±0.009 


4. 


c 


->* + ... 


57.8 


0.541 ±0.036 


5. 


c 


->* + ... 


91 


0.522 ± 0.022 


6. 


b 


->*+... 


29-35 


0.789 ±0.022 


7. 


b 


->*+... 


91 


0.699 ±0.009 



Table 1: Experimental measurements of (xq) (W,M) 

Errors have been evaluated by taking statistical and systematic errors in quadrature. First 3 
entries stand for the direct production of D* mesons at different centre-of-mass energies; the last 
4 data for the mean quark energy have been extracted by unfolding the inclusive lepton (e, fj) 
spectra from heavy Q decays. 

As mentioned above, the PT approach advocated in this paper can not pretend to fully 
describe exclusive heavy hadron spectra (with D* an example). Our treatment of the hadroniza- 
tion stage that implicitly appeals to duality arguments makes it plausible to rather apply this 
approach to inclusive quantities such as the lepton energy distributions. Nonetheless, for lack 
of anything better, we take the measured mean energy of D* as a representative of (xq) to be 
compared directly with the PT motivated prediction for the quark energy losses. 

The preliminary analysis has shown fffl that an independent fitting of the D and L(epton) 
data results in the best-fit curves Ao(A) and At, (A) that cross just at the best-fit value of A. 
This was the observation that motivated us to look upon af as a process independent quantity 
to confront the c and b measurements 1-7 with a unique one-parameter PT prediction^. 

Hereafter we fix heavy quark masses to be 

M c = 1.5 GeV , M b = 4.75 GeV . (4.1) 
(Sensitivity to the 6-quark mass will be discussed below.) 

4.1 Fitting mean quark energies. 

Fig.H shows the quality of the fit to 7 separate data of the Table [L] together with the total x 2 as 
a function of A^ 5 ^ in the G^-model with C2 = 2.86 (the best-fit value). Some explanation is in 
order. In this Figure (and similar plots below) for each datum the ratio is displayed 

theor. - exp. 
exp.error 

9 in spite of the fact that such a naive approach suggests the same theoretical expectation for the two physically 
different data # 3 and # 5. 
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against the right vertical scale. Dashed horizontal lines mark la levels for each single datum. 
The two curves (which should be read out against the left scale) show the squared deviation of 
the points (|4.2|) from the median, that is, total x 2 ■ The solid curve sums up all 7 data, while the 
dash-dotted one accumulates the "high W" data only. The "high W" sample we define excluding 
the entries #1 and #6 which correspond to the quark mass-to-energy ratios 

M/E « M c /5GeV ~ M 6 /15GeV « 1/3 . 

These two entries are subject to significant non-relativistic corrections. The fact that the two fits 
are consistent is good news: it shows that the non-relativistic effects have been properly taken 
into account^ in the PT radiator ( j2.7| ) . 

A 0.191 




0.45 0.5 0.55 0.6 0.65 0.7 

A (5) (GeV) 

Figure 2: A-dependence of the F-model fit to mean energy losses (two-loop af with n/ = 5). 
The right scale shows normalized deviation between a theoretical prediction and an experimental 
datum ( [4.21) . Dashed lines mark the la band. Solid and dash-dotted lines show the values of x 2 
(against the left scale) for all data and the high-VT data sample correspondingly. 

Fig.|3| demonstrates sensitivity of PT description to the "confinement parameter". Here we 
have chosen the G2~ m odel for a change. The first thing to be noticed is that with the best-fit 

10 the relativistic version of (2.4-7) reported earlier failed to properly embody the "6 at 32" datum # 6 
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parameter C% = 2.86 one obtains the same value A^ 5 ) « 580 MeV, comparable quality of the fit 
Xmin ~ 0.7 and even the same dynamics of each of 7 data as in the above F-model description 
(Fig.§). In the upper part of this plot two marginal values of C<i are also shown which correspond 
to one standard deviation from the total 7— fit: x^ in (2.53) = x^ in (3.21) = x^ in (2.86) + 1. 



X B(91) O B(32) □ C(91) + C(58) « D(9l) * D(30) X D(10) 



X All data 



X High W data only 



x X X x x 




0.4 0.6 
l (5) 





"1 1 




A (5) = 0.50 : 




\ A (5) =0.58 




f\ A (5) =0.66 ; 




ul 1 , - 


0.5 


1 5 10 




k (GeV) 



A w (GeV) 

Figure 3: On sensitivity of the G 2 -model fit to the shape of a"f in the origin (two-loop, n/ = 5). 
The values of C% in the upper plots correspond to one standard deviation from the best-fit 
(bottom-left). Bottom- right graph displays the margin in low momentum behaviour of a e f. 



In Figj^ comparison is made of the quality of the fits within different models for the effective 
coupling (two-loop, nj = 5). For each value of A the A/C parameter has been adjusted to 
minimize the error (one-parameter fit). The upper scale shows corresponding values of at 
LEP recalculated from a e f with use of the relation (2.17b). 
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It is worthwhile to notice some peculiarity of the Gi-model. This model is "too soft" in a 
sense that it induces the negative preasymptotic power term oc k~ 2 in af(fc), which correction 
suppresses af in a relatively high momentum region. 

Fig.|| illustrates this peculiarity. Here the couplings corresponding to the best-fit A/C values 
for A^ 5 ) = 580 MeVare compared. Solid lines (F, G2-4) correspond to x 2 ~ 0.7. In the Gy- 
model shown by dash-dotted line (x 2 ~ 1.2) «f stays noticeably smaller above 1.5GeV before 
the perturbative logarithmic regime sets up and all the models merge. 



0.120 



«-(M 7 ) 
0.125 MS V z/ 0.130 



«f(k) 



0.135 




0.5 



0.0 



Two-loop a s etf (n f =5) 



0.4 



0.5 0.6 
A (5) (GeV) 



0.7 




k (GeV) 
it of fi 

model underestimates radiation above 1.5GeV. 



Figure 4: One-parameter fits to energy Figure 5: Best-fit af for A^ 5 ) = 580MeV. G x - 
losses within different models for a e f. 



As a result, to compensate for reduced radiation intensity the best-fit A value (and thus 
a(M z )) within the Gi-model tends to be larger compared to "sharp" models F, G 2 , ■ ■ ■ 

Leaving G\ aside we conclude that both the quality of the fit and the value of A the "sharp" 
models point at, hardly exhibit any model dependence. From Fig.[| (see also Fig.|3|) we deduce 

A (5) 



580 ± 80 MeV . 
Being translated into the MS parameter this gives 

«ms(M k ) = 0.127 ±0.003. 



(4.3a) 



(4.3b) 



The error here is purely statistical (one standard deviation). 

In what follows we shall be using the 0.003 shift in aj^(M z ) induced by the Gi-model (see 
Fig.f| and Table 2 below) as a rough estimate of systematic uncertainty due to possible "soft" 
preasymptotic power effects in the running coupling. 
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Fig.|| demonstrates consistency of the total G2~model fit with fits to various subsets of data: 
D* (items 1-3 of Table 1), "leptons" (4-7), high W (items 2-5, 7). Here x 2 is plotted against 
the reference value a^g(M^). 

D* data are more restrictive since they have smaller experimental errors than inclusive lepton 
measurements. Low-jy points b — depton+. . . at 32 (#6) and, especially, D* at 10 (#1) are quite 
important as they provide lever arm for scaling violation. Inclusion of these two measurements 
does not spoil the fit, x 2 /5 (total fit) ~ X 2 /3 (high-W 7 ), but increases its quality reducing 
statistical error by factor 2. 



3.0 




0.0 



0.12 



0.125 



0.13 



0.135 



a ^( M z) 



Figure 6: On consistency between the fits to different subsets of data on mean energy losses 
(two- loop, Uf — 5). Horizontal lines mark one standard deviation levels. 



4.2 Energy ratios and A determination. 

In Fig.0 results of the G^-rnodel fit to the ratios of (xq) at different energies are shown. 
Solid curves accumulate the squared errors in the description of four ratios: D*(91)/D*(10), 
D*(91)/£>*(30) and two ratios from leptonic quark decays, C(91)/C(58) and 5(91)75(32) As 
we have discussed above, one expects such ratios to be protected against our ignorance about 
the confinement physics. Indeed, a rather high stability in the quality of the fit inside a huge 
range of variation of the "confinement parameter" is seen. Bottom-right insertion displays 
corresponding shapes of af. Within the chosen interval of C2 the characteristic value of A (|4.6|) 
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varies from 0.09 (C 2 = 10, A = 0.7) up to 0.41 (C 2 = 1.4, A = 0.5) that is changes by factor 2 in 
both directions around the value 0.19 ( f4.6b|) that we have obtained describing absolute energy- 
losses. 




Figure 7: (^-dependence of energy ratios (Cr2-model, two-loop, rif = 5) with x 2 for the first four 
ratios shown by solid lines and correponding shapes of of (bottom-right). The ratio 5(91) /C(91) 
not belonging to the fit exhibits strong (^-dependence, contrary to generic c/c and b/b ratios. 

Also shown in Fig.[7] is the mixed bottom-to-charm ratio £?(91)/C(91). This one does not 
belong here and, contrary to the generic c/c and b/b quark ratios, strongly depends on C2 as 
expected. 

We observe that three of four generic ratios show up no variation with C2 at all. It is the ratio 
D*(91) I D*(10) only that induces some negative correlation between A and A: the latter moves to 
larger values with decrease of low-scale interaction intensity. Such a systematic drift is natural: 
charm quark mass is too small to completely protect the normalization point D*(10) against 
W - dependent (sic!) confinement efects at total energy as low as 10 GeV. At the same time the 
first ratio dominates in the fit while experimental accuracy of three others is not sufficient at 
present to provide a direct safe way of measuring the A parameter. 

Roughly one might present the result of fitting quark energy ratios as 

A (5) = 0.60 ± 0.15 (stat.) ± 0.10 (syst.) . (4.4a) 
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If statistical and systematic errors in Fig.ff] were uncorrelated (which is not the case), (|4.4a|) 
would correspond to 



a m (M z ) = 0.128 ± 0.007. (4.4b) 

For the time being it will suffice to conclude that determination of A*- 5 -* from quark energy ratios 
is consistent with that from absolute energy losses, cf. ( [4.3| ). 

4.3 On bottom quark mass. 

Fig|| illustrates sensitivity of the total 7-fit to the bottom quark mass. Here we have fixed 
M c = 1.5 GeV and looked for minimal x 2 with respect to variation of A/C2 (1-parameter fit) for 
given A and M&. 




0.4 0.5 0.6 0.7 

A (5) (GeV) 

Figure 8: Mf, dependence. Fits within F and G2 models are shown (two-loop, nf — 5.) 

Quark mass dependence of the spectrum (|2.4|) - (|2~77D and, thus, of energy losses is basically 
logarithmic (apart from non-relativistic corrections (M/W) 2 and "confinement" effects [i a /M). 
Nevertheless, a certain range of "preferable" bottom quark masses may be read out from Fig.|^. 
Within one standard deviation (G 2 -model) 

M b = 4.73 ± 0.38 GeV. (4.5) 

A better understanding of the nature of the mass parameter M as it enters PT formulae is needed 
before such an analysis may be used for "measuring" bottom quark mass. 
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As for now, one finds satisfaction in noticing that the value M& = 4.75GeV we have been 
using throughout this paper does fit nicely into the PT picture of mean energy losses. 

4.4 Integrated coupling as invariant of the fit. 

Being free to play around with the detailed shape of a e f in the origin, one finds it necessary 
at the same time to fix some characteristic measure of the radiation intensity in the non-PT 
momentum region. Fig.|5] suggests trying an area under the curve for an invariant parameter of 
the fit. The result is quite impressive: areas under the G2-4 an d F couplings are practically 
indistinguishable. Introducing 

A(u) = - dk^Q, (4.6a) 
/1 Jo IX 

we get (with one standard deviation error estimated from the G2~niodel margin, see Fig.^J) 

A(2GeV) = 0.190 ±0.010. (4.6b) 

This integral measure not only proves to be quite stable against the choice of /ow-momentum 
regularization but also resistant to different approximations for the high- momentum tail of a e f(k). 

Table 2 accumulates charateristics of various 7-fits and demonstrates an amusing stability of 
the value ([4.6b]) . This justifies the qualitative expectation of subsection 3.1 that it is the integral 
of the coupling as a characteristic measure of confinement (hadronization) effects in inclusive 
energy spectra, \i a of eq.(|3.9|), that is responsible for the low-momentum contribution to mean 
energy losses. 

Thus we find empirically that the characteristic integral ( |4.6[ ) turns out to be a fit-invariant 
quantity which one has to keep fixed to describe the absolute values of energy losses. As pointed 
out by V.N. Gribov, it can be looked upon as the long-distance contribution to the QCD field 
energy of a heavy quark. It is worthwhile to notice that such an integral appears in the relation 
between the running heavy quark mass at scale // and the pole mass @ 

MP0le_M W = ^/ ^L'^ = C F rdn'^-^C Flla . (4.7) 

3 J\k\<n {ZTT) 6 K z Jo 7T 

4.5 Two-loop a e s a with heavy quark thresholds and a-^r determination. 

It is important to notice that the different approximations for the high-momentum tail of a e f 
listed in Table ^provide similar quality fits and preserve the A value but, at the same time, lead 
to systematically different values of ay^(Mz)- Fig.|9] helps to relate the values of A parameter 
for different approximations for the running a e f to the reference value oe^g(Mz)- 

In the problem under consideration one probes the coupling at Mz scale only indirectly. 
For the first thing, "half" of the data belong to smaller Ws. Moreover, even when the LEP 
measurements are concerned, the main contribution to energy losses originates from a broad 
logarithmic integral of a c f(k) running from k < Mq up to k < Mz, so that momenta just above 
Mft (as well as between M c and M&) play quite an essential role. The reference values of a^g 



24 



Model for 


Best-fit 






y 2 

A 


Integral A 


of 


parameter 


A (MeV) 


OiTIaiM?) 


5 d.o.f. 


( 


4.6b 


) 


2-loop, 


A— 0.184 


730 ± 95 


0.125 


0.66 


0.184 


with 


C 2 = 2.48 


725 ± 85 


0.125 


0.64 


0.191 


c, b 

1 


C 3 = 4.41 


710 ± 80 


0.124 


0.66 


0.193 


thresholds; 


C 4 = 7.79 


705 ± 80 


0.124 


0.68 


0.194 


A( 3 + 2 ) 


Ci= 1.38 


820 ±100 


0.128 


0.69 


0.189 


2-loop, 


A = 0.190 


585 ± 85 


0.127 


0.73 


0.187 


with 5 


C 2 = 2.88 


585 ± 80 


0.127 


0.68 


0.191 


massless 


C 3 = 5.59 


575 ± 80 


0.127 


0.69 


0.193 


quarks; 


C 4 = 10.8 


575 ± 75 


0.127 


0.70 


0.193 


A (5) 


Ci= 1.46 


655 ± 95 


0.130 


0.75 


0.190 


1-loop, 


A= 0.193 


480 ± 70 


0.120 


0.71 


0.18£ 


\ 


with 3 


C 2 = 2.22 


485 ± 65 


0.120 


0.66 


0.190 


massless 


C 3 = 3.94 


480 ± 60 


0.120 


0.66 


0.191 


quarks; 


C 4 = 6.90 


480 ± 60 


0.120 


0.66 


0.191 


A (3) 

1-loop 


Ci= 1.20 


515 ± 75 


0.122 


0.69 


0.189 


1-loop, 


A = 0.209 


310 ± 50 


0.133 


0.90 


0.189 


with 5 


C 2 = 2.68 


315 ± 50 


0.133 


0.79 


0.191 


massless 


C 3 = 5.48 


315 ± 50 


0.133 


0.78 


0.191 


quarks; 


C 4 = 10.9 


315 ± 50 


0.133 


0.79 


0.191 


A (5) 

A l-loop 


Ci= 1.28 


330 ± 50 


0.134 


0.84 


0.189 



Table 2: Best 7— fits within various models for effective coupling. A/C, a^g, x an d A are given 
for the central A values. 

appearing in Table |2| emerge as a result of extrapolation from intermediate momentum scales that 
dominate in the fit. Such an extrapolation is sensitive to details of high- momentum behaviour 
of the running coupling. When n/ is taken smaller and/or the two-loop effects are being taken 
into account, a e f becomes a steeper falling function of momentum and the resulting value of 
ajfis(M z ) decreases. 

Even an account of heavy quark thresholds (which makes a c s s (k) a steeper function below 



k < Mb) drives down the ajjc value. To demonstrate this effect we include Figs. 10 and 11 



showing quality of the PT description of absolute energy losses with account of the second loop 



and c, b quark threshold effects in the running coupling. In Fig.|i0| model dependence of the 
total 7— fit is shown (cf. Fig.£|); Fig.|Tl| collects fits to different subsets of data (cf. FigH). 

Our conclusions about relative stability against the choice of low-momentum regularization 
model and about consistent description of leptonic, D* and high-PF data hold. From these plots 
one obtains (see also Table 0) 

A (3+2) = 72Q ± ogo ( gtat ) ± Q15 ( sygt ) GeV ^ ( 4ga ) 
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Figure 9: Relation between A and oe^g(Mz)' one-loop (dotted), two-loop (dashed) for 3 and 5 
massless quarks and the two-loop coupling with c, b quark thresholds (solid). 



where we have singled out systematic uncertainty due to model dependence. This results in 

a m (M z ) = 0.125 ± 0.003. (4.8b) 

It is important to stress that the second loop effects are there in the running coupling and rif is 
not a free parameter: quarks (and their masses) are what they are. Therefore a huge interval of 
a MS vames seen m Table ^| has nothing to do with actual theoretical uncertainty in determining 
this important datum. 

Before we turn to discussion of systematic errors some comment is in order. It has to do with 
the number of acting quark flavours and will eventually give us an additional consistency check 
of the approach. 

An attentive reader could have noticed that using 5-flavour running coupling seems to con- 
tradict the very logic of the present paper. Indeed, in Section we constructed the energy 
distribution corresponding to direct production of a heavy quark Q. To do so we disregarded sea 
mechanism of heavy quark production and omitted the second loop anomalous dimension term 
related to QQQQ final states. 

However, there is yet another subtle source of copious QQ pairs, namely, the running coupling 
«f embodied into PT radiator. The way the coupling runs in theoretical expressions for inclusive 
characteristics depends on experimental setup: Veto on fermion pair production in the final state 
suppresses effective interaction strength at large momentum scales (QCD coupling decreases 
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A (3+2) (GeV) 

Figure 10: Quality of the total fits versus 
A (two-loop af with c, b thresholds). 
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Figure 11: On consistency of G^-model fits to 
subsets of data (two-loop; c, b thresholds). 



faster; an increase of a em slows down). Therefore, to preserve the logic of the approach, we 
better make it clear that the use of nf = 3 effective coupling which goes along with suppression 
of additional cc, bb pairs is consistent with the reported result ( f4.8b| ) . 



In Fig.|T^ G2-model couplings are shown with the best-fit parameters A and C2 listed in 
Table || As far as low momenta are concerned, variations due to n/ and one vs. two loops 
should be looked upon simply as different models for trial regularized coupling. All of them 
do the job. In particular, the 3-flavour af which interests us at the moment (the last line of 
Table |3|) does provide an excellent fit, x 2 /d.o.f.~ 0.6/5. 



G^-rnodel 


"/ 




c 2 


x 2 


A(2GeV) 


af(2GeV) 


af(M b ) 


afiM z ) 




+thresholds 


3+2 


725 


2.48 


0.65 


0.191 


0.530 


0.303 


0.135 


0.125 




5 


315 


2.68 


0.79 


0.191 


0.443 


0.302 


0.145 


0.133 


1-loop 


4 


395 


2.43 


0.72 


0.191 


0.463 


0.303 


0.139 


0.126 




3 


485 


2.22 


0.66 


0.190 


0.488 


0.306 


0.134 


0.120 




5 


585 


2.88 


0.69 


0.191 


0.494 


0.303 


0.138 


0.128 


2-loop 


4 


770 


2.72 


0.64 


0.190 


0.549 


0.305 


0.131 


0.120 




3 


935 


2.56 


0.62 


0.190 


0.608 


0.305 


0.124 


0.112 



Table 3: Characteristics of the best G2~ m odel fits. 
High stability of the A value (4.6b) is confirmed once again. Moreover, one more interesting 
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0.5 1 5 10 



k (GeV) 

Figure 12: Best-fit G2~iriodel effective couplings, rif = 5,4,3 progress from left to right. 



fit-invariant quantity emerges, namely, the value of the effective coupling at the bottom mass 
scale0, 

af(M b ) « 0.30. (4.9) 



The curves in Fig.[L2| are quite different below 2GeV, focus around 5GeVand diverge again at high 
momenta. The columns adjacent to the double-lined one in Table § illustrate this behaviour. In 
particular, the 3-flavour coupling continued to the LEP scale is substantially lower compared to 
the "3+2" coupling we have been using before (0.124 vs 0.135). However, it would be erroneous 
to read out Oj^g = 0.112 from the last line of Table |^. 

As we see, heavy flavours in af are indeed irrelevant for the problem under consideration. 
However, they do contribute to the evolution of the standard QCD coupling. Therefore, having 
found the magnitude of a e f at intermediate scales k ~ 5GeV which dominate in the problem, one 
has to apply the 5-flavour evolution from ~ 5 up to 91 GeV aiming at extraction of the reference 
value ajj^(Mz). 

Since the starting value a^ ff (M{,) = 0.305 practically coincides with that from the previous 
"3+2" analysis (0.303), so do the results. Choosing the starting point of high momentum ex- 



ii 



since magnitude of a e s s around 5GeV S> A is insensitive to A/C and depends only on A, the same conclusion 



(4.9) follows from F and G3.4 models 
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trapolation in between k = Mb/2 and k = 2Mb leads to 



MS 



0.125 ± 0.003. 



(4.10) 



The, slightly inflated, error here is moderate, and the central value perfectly matches with the 
result of the previous analysis (|4.8b|) based on "3+2" coupling. 



4.6 Systematic errors and prospects. 

First of all, there is a problem of correspondence between theoretical prediction and the data. 
As we have stressed above, formulae of this paper were designed to describe mean energy of a 
primary quark produced in the e + e~ annihilation vertex. We were treating experimental numbers 
as corresponding to direct QQ production. We do not feel in a position to judge to what extent 
the experiments actually met such an expectation. 

The good news is however that the bottom sector is safe in this respect: g — > b sea component 
is vahishingly small at present, and so is charm production at low energies il. Therefore only 
LEP charm data should concern us here. Bearing this is mind it is important to notice that the 
D* LEP datum ^3 (which is the most precise measurement and therefore practically the only 
one vulnerable) does correspond to primary charmpj MC modeling was used to subtract the 
gluon component (with an uncertainty included into experimental systematic error). Reliability 
of such a subtraction has been recently verified by the first OPAL measurement of the charm 
production via gluons @. 

Apart from this, one can think of the following sources of systematic error of a-^(Mz) 
determination. 

higher orders. Estimate: 0.002; Comment: optimistic. 
Follows from analysis of approximate relation ( |2.17b| ) between physical coupling af and a^fg 



and of forceful exponentiation of a s terms in the radiator (|2.7|) . Apart from three-loop analysis, 
an exact O (a 2 ) theoretical calculation of (xq) would be helpful. 

power effects in a e f(k > 2GeV). 0.003; arbitrary (hopefully conservative). 
Is based on comparison with the "soft" Gi-model (see above). Variation within "hard" regu- 
larization models (F, G2-4) is below 0.001 (see ( [4.8|) ). To gain quantitative theoretical control 
over the 1/k 2 power term in the effective coupling and, thus, to reduce corresponding systematic 
error, might be not as hopeless a goal as it seems. 

kinematical effects. 0.000; safe. 
The structure of the PT spectrum is such that it respects the kinematical boundary x > x m i n = 
2M/W only in the first a s order. Contribution to (xq) from potentially dangerous region between 
Xm\n and 2xmiTi 

can be estimated as 



A {x Q ) ~ a s (W)/n ■ (1 - (x Q )) ■ x 2 min , 



which value is well below 1% even for x min as big as 1/3. 

12 An older average = 0.494 ± 0.014 which we have used in the previous analysis §1 was contaminated 

by the sea. We are indebted to P.Mattig for clarifying this point 
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quark masses. 0.002; conservative. 
Given that the structure of low-momentum contribution proportional to A naturally embodied 
into PT analysis reminds that of the shift from Euclidean to "on-shell" quark mass (see (|4.7|) ), 
one may worry about double counting if the pole mass is used for M in the PT formulae. The 
problem should be studied theoretically. Meanwhile, let us mention that one may get equally 
good description of absolute values of mean quark energy losses with smaller (Euclidean?) quark 
masses plugged in. For example, by taking 



M c = 1.30 GeV; 



M h = 4.50GeV 



one obtains^ 



G2-model 


c 2 


x 2 


A(2GeV) 


a Ms( M z) 


A^ =600 


3.60 


0.72 


0.166 


0.1280 


A (3+2) =74Q 


2.94 


0.66 


0.166 


0.1257 



We observe that the resulting value a^(Mz) hardly increases by 1 per mil. 

Another feature of the alternative set of quark masses worth noticing is a systematic decrease 
of the value of characteristic integral A. This may be a welcome trend bearing in mind a recent 
analysis of power corrections to jet shape variables, see [35]. For the time being we choose to 
look upon this 10% shift as a systematic uncertainty which will be greatly reduced when proper 
theoretical understanding of the nature of the quark mass parameter is achieved. 

Having said that we modify ( |4.6b|) and present the final result for the integrated coupling as 



A(2GeV) = 0.18 ± 0.01 (stat) ± 0.02 (syst) . 



(4.11) 



Finally, in Fig.|T^ we demonstrate the W^-evolution of realistic c- and b-quark spectra obtained 
within the F-model with one-loop 3-flavour af. These curves correspond to the values of A and 
A providing the best common fit to mean energy losses as described above. For comparison the 
best-fit G^-model curves are also shown for LEP energy. F- and G-model spectra are quite close 
to each other. We may conclude that is suffices to fix the integral parameter ( [4.6 b[ ) together with 
the value of A to predict the differential energy distributions with a reasonable accuracy! 1 ^. 

A stable numerical procedure for numerical evaluation of the inverse Mellin transform 
remains to be designed. 



5 Conclusions 

Jets initiated by heavy quarks (b,c) are now extensively studied experimentally. The interest 
to this subject is connected not only with testing the fundamental aspects of QCD but also 
with its large potential importance for measurements of heavy particle properties: lifetimes, 
spatial oscillations of flavour, searching for CP violating effects in their decays etc. Properties 



13 4.5 is a preferable bottom mass value for M c = 1.3, analogously to Mb — 4.75 for M c =1.5; see abo ve, F ig ]8| 
14 Let us notice that the differential quark energy spectra obtained by the inverse Mellin transform (2.4), ( |2.6[ ) 
violate kinematical boundary D{x<2M/W) =0 at O (oq) level 
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of b-initiated jets are of primary importance for analysis of the final state structure in the ti 
production processes. 

In this paper we presented results of the study of the inclusive energy spectra of leading 
heavy-flavoured particles (Hq) based on the perturbative expression ( |2.4jJ - ( |2.7| ) for heavy Q 
distributions that emerge after taking into proper account multiple gluon bremsstrahlung off 
the QQ pair. 

Our approximation includes the two-loop anomalous dimension, keeps track of the collinear 
logarithms a In W and a In M, soft double- logarithmic a ln 2 (l— xq) and essential single-logarithmic 
aln(l— xq) contributions in all orders. At the same time, it embodies the exact first order 
result 0(a) for the inclusive energy distribution which property is necessary to account for non- 
relativistic suppression of gluon radiation in the situation when jet energies are comparable with 
a heavy quark mass. 

The spectra we have been considering correspond to direct Q production. Such an approach 
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disregards the sea contribution to heavy quark yield which is numerically small at present and 
foreseeable energies. Accordingly, a specific term in the non-singlet two-loop anomalous dimen- 
sion has been omitted which is due to QQQQ final states and would violate the "multiplicity" 
sum rule ( |2.1|) . 

Apart from this simplification, the formal relative accuracy of the perturbative result ( |2.4| ) - 
( p.7|) is O (a 2 (W 2 )) which estimate is uniform in xq. Perturbative corrections at the quark mass 
scale which were left out of control start from O (a 3 (M 2 )). 

We have formulated the perturbative result ( |2.4D - ( ^.7[ ) in terms of a "physical" coupling, 
the one which directly measures radiation intensity of relatively soft gluons O. This coupling 
constant is different from (roughly, 10% larger than) cty^ and is approximately related to the 
latter by (|2.17b| ). Expressed in terms of af, the radiator ( |2.7| ) acquires the two-loop contribution 



A^ 2 ) which is free from an artificial n^-dependence i! and proves to be numerically negligible. 

There are two ingredients of the standard approach to description of Hq spectra. Here one 
starts from the phenomenological non-PT fragmentation function for the Q — > Hq transition 
and then traces its evolution with increase of the annihilation energy W by means of PT QCD. 
Being formally well justified for describing the IV-evolution, this approach however leaves the 
dependence on the quark mass, M, basically out of the PT control.^ 

Motivated by the LPHD concept, we attempted to derive pure PT predictions without in- 
voking the phenomenological fragmentation function. Starting point for such an attempt was 
the observation that an appearance of the parton model motivated peak in the non-PT frag- 
mentation at large xq 00 can be attributed to the Sudakov suppression effects provided one 
feels courageous enough to continue the PT description down to the region of gluon transverse 
momenta, k± ~ A • (W/ V MA) 0,2 , which at present (and foreseeable) energies looks dangerously 
close to the non-perturbative domain. 

When getting rid of the transverse momentum cutoff one faces the problem of the formal 
"infrared pole" in a s and is forced to introduce the effective non-singular coupling ajf(/c) that 
remains finite at k<A. It is not easy to justify the very notion of a s at small momenta where the 
PT quark-gluon language seems to be hardly applicable at all. In the problem under consideration 
it can be related to the effective measure of intensity of accompanying particle production at the 
confinement stage of the Hq formation, a finite number of light hadrons produced in addition to 
the (W dependent) particle yield due to PT-controlled gluon bremsstrahlung at the first stage 
of the hard QQ pair creation process. 

We have checked that, in accordance with expectations based on general factorization prop- 
erties, our ignorance about confinement does not affect the W^-dependence of the mean energy 
losses (xq) (W, M). Thus the ratio (xq) (W)/ (xq) (W ) can be used to extract the fundamental 
QCD parameter A. Another infrared stable quantity found empirically is the PT prediction for 
the normalized peak position, XQ ak (W)/x 1 Q ak (W ). 

At the same time, the absolute values of (xq), peak positions and particle distributions in 
general show up substantial variations with the value (momentum dependence) of af in the 
origin chosen as an input for calculating the PT radiator ( |2.7|) . We have compared experimental 
values of (xq) extracted from leptonic c- and b-quark decays and (xd*) directly with the PT 



15 A plausible rescaling procedure for extracting the non-PT part of the B meson fragmentation from the existing 
data on D meson spectra has been suggested recently in 
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motivated prediction for the quark energy losses. 

Our treatment of the hadronization stage that implicitly appeals to duality arguments makes 
it plausible to apply purely perturbative approach to inclusive quantities such as the lepton 
energy distributions. Meanwhile, it can not pretend to describe exclusive heavy hadron spectra 
such as those of D*, for example. Therefore, taking the measured mean energies of D* for (xq) 
might introduce some systematic uncertainty to the fit, the error which is difficult to estimate. 

Nonetheless, such a comparison has demonstrated that radiative charm and beauty losses 
may be consistently described in the energy range W = 10 — > 90GeV within a variety of models 
for low momentum behaviour of the infrared finite effective coupling af. By tuning the shape 
parameter of a model one achieves an accuracy of the common fit as good as x 2 /d.o.f.~ 0.7/5 (7 
data points; two free parameters, A and A/C). 

We have checked that the fit to all data is consistent with the fits to high-W^ , D* and inclusive 
lepton data samples. 

From the first sight, we have gained not much profit since we had to pay for eliminating 
an arbitrary non-PT fragmentation function by introducing another arbitrary function a e f(k) at 
k < 2GeV. An essential physically important difference between two approaches is that a c f as a 
key ingredient of the PT-LPHD approach is supposed to be universal so that (quite substantial) 
differences between inclusive spectra of c- and b- flavoured hadrons should be under control. By 
simultaneous fitting of the available experimental data on the energy losses in charm and beauty 
sectors we have demonstrated consistency of the hypothesis of a e f universality. 

Moreover, the same notion of the infrared-finite effective coupling can be tried for a good many 
interesting problems in the light quark sector. An incomplete list of such phenomena for which 
the PT analysis has been carried out recently to next-to-leading order includes transverse @ 
and longitudinal momentum distributions in hadron-initiated processes, the energy-energy 
correlation §1, the thrust @ and the heavy jet mass distribution H in e + e~ annihilation. 

The presence of the exponential of the characteristic integral over gluon momenta which 
emerges after all-order resummation of the Sudakov logarithms 0E§ is a common feature of 
the corresponding PT expressions. It is straightforward to derive quantitative PT-motivated 
predictions by implementing the universal a e f in these integrals and, at the same time, by getting 
rid of non-PT "hadronization" effects which are usually taken into account by convoluting the 
PT distributions with phenomenological fragmentation functions (initial parton distributions for 
the case of hadron-initiated processes). 

Carrying out this laborious but promising program one should get a valuable information 
about the confinement physics as seen through the eyes of the integrated influence of the large- 
distance hadron production upon inclusive particle distributions and/or event characteristics. 

To this end the numerical results of our analysis of the heavy quark energy losses, namely, 

i f 2GeV a ett (k) 

(2GeV)- 1 dk sV ' = 0.18 ± 0.01 (exp) ± 0.02 (theor) , (5.1a) 

JO IX 

a m {M z ) = 0.125 ± 0.003 (exp) ± 0.004 (theor) , (5.1b) 

should be looked upon as a first hint for the more detailed study including the light quark 
phenomena discussed above. 
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One comment is in order concerning the nature of the results (|5.1[) . The parallel description 
of c,b losses is sensitive to a e f in the low momentum range, say, of the order of (and below) Mb, 
so that the integral characteristic ( |5.1aJ) gets fixed quite sharply (theoretical error is due to our 
ignorance about the heavy quark masses and should be eliminated in the future). 

At the same time such a description proves to be quite liberal to details of the high-momentum 
behaviour of the coupling (one- versus two-loop a s , number of active quark flavours). From this 
point of view the value (|5.1b|) that one extracts in the end of the day is rather a tribute to 
(unfortunate) tradition. a^(Mz) should be looked upon as a reference value which emerges 
after theoretical extrapolation rather than a quantity that is "measured" by the above analysis. 
Theoretical error in ( |5.1b|) is dominated by potential k~ 2 power correction effects in a e f(k). 

Let us stress again that in the low momentum region behaviour of the effective coupling is 
poorly known not because of the limited knowledge of higher order effects but rather because 
of an essentially different physical phenomenon that enters the game, the one that is usually 
referred to as confinement. 

From this point of view our notion of the infrared finite af differs from one that emerges from 
the three-loop analysis based upon the renormalization-scheme-invariant approaches to the e + e~ 
annihilation cross section @@ and the r-lepton hadronic width £3. Nevertheless it is worth 
mentioning that the value of the couplant — (0) and the integral measure obtained in [45| in the 
framework of the Minimal Sensitivity Principle il are consistent with (|5.1a|). 



Phenomenological verification of the fact that the effective QCD coupling stays numerically 
small would be of large practical value. Gribov theory of confinement @ demonstrates how 
colour confinement can be achieved in a field theory of light fermions interacting with compar- 
atively small effective coupling — a fact of potentially great impact for enlarging the domain of 
applicability of perturbative ideology to the physics of hadrons and their interactions. 
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Appendix A : Running coupling in the Radiator 



A.l Notation and Kinematics 

We consider production of a QQ pair with on-mass-shell 4-momenta pi, p 2 and a gluon with 
momentum by a colourless current with the total momentum q^. Hereafter for a sake of 
simplicity we measure all the momenta and the quark mass in units of the total annihilation 
energy (q 2 = W 2 = 1). 

In terms of quark and gluon energy fractions, 

Zi = 2(piq) , z = 2(kq) ; z 1 + z 2 + z = 2 ; 

2p!p 2 = (q - k) 2 - 2m 2 = 1 - z - 2m 2 + k 2 = z\ + z 2 - 1 - 2m 2 + k 2 ; 

(A.l) 

virtual quark propagators are 

«x = (px + k) 2 - m 2 = 1 - z 2 , k 2 = (P2 + fc) 2 - m 2 = 1 - z x . (A.2) 

(Relations dAT2|) do not imply k 2 = 0.) 

Three-body kinematics restricts the difference of quark energy fractions as 

A\k\ 2 -P 2 > ( Zl -z 2 )\ (A.3a) 

where (3 = )3(z) represents the quark velocity in the rest frame of the QQ pair (QQ cms): 

^9 Pali 4m 2 Am 2 

(3 2 ^^ = 1- = 1 - — . A.3b 

E qq {P1+P2) 2 (q-k) 2 

A. 1.1 k 2 = 

For the on-mass-shell gluon case one has 4\k\ 2 = z 2 ; (q — k) 2 = 1 — z and (|A.3|) gives 



1 



z\ - z 2 



(A.4) 



z, G c Basis. When the gluon energy is kept fixed, it is convenient to use the scaled difference 
of quark energies as a complementary variable related to the gluon angle in the QQ cms: 

_ Z1-Z2 2(p l ~p 2 )k 

u= = — ■ — = /3cos6 c ; -p < u < +(3. (A. 5) 

z 2(p 1 +p 2 )k 

The maximal value \u\ =/3 is reached when the gluon is collinear with one of the quarks. 
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Introducing maximal quark velocity v, 

/3 2 <v 2 = I- Am 2 > z , (A.6) 
one may present the integration phase space as 

rv 2 r /3 



J dz\ J dz 2 — | J zdz J du . (A. 7) 



Useful relations: 



l-u 2 = 4(1 - zi)(l - z 2 ) ■ z~ 2 ; (A.8a) 
1 — z-t 1 — Zo\ z 2 4 , , , . 

+ 1 ' = 71 SH V - 2 = ~, 9 ~ 2 ■ ( A - 8t) ) 



1-Z 2 l-zj {1- Zl ){l-Z 2 ) ' l-U 2 

x, y Basis. For the purpose of deriving single-inclusive quark spectrum we break the symmetry 
between zi, z 2 and denote 

x = z x \ y = l-z 2 . 

One has to fix x and integrate over y in the limits y_ < y < y + which follow from the kinematical 
restriction ( |A.4| ). In terms of x, y the latter takes the form 



y± = — r- (x — 2m 2 ± V x 2 — 4m 2 ) . 



y 2 {l -x + m 2 )- y(x - 2m 2 )(l - x) + m 2 (l - x) 2 < . (A.9) 

This gives 

" 2(1 -x + m 2 ) 
Introducing 

z Q = - (x - 2m 2 + Vx 2 - 4m 2 ) = x + O (m 2 ) , (A.10) 
one may write (cf. (|2.8| )) 



y_ = -m 2 = -m 2 ■ (l + O (m 2 )) 



(A.ll) 



z x 

Useful relations: 

1 — x 



y+-y- = a/x 2 - 4m 2 , (A. 12a) 

1 — x + m z 

(y-)- 1 - {y + Y l = Vx * " 4m2 m~ 2 , (A.12b) 

1 — X 

(y + ) 2 -(y^) 2 = (x-2m 2 )Vx 2 -4m 2 ( _ i l ~ X V. (A.12c) 

V 1 — x + m z J 
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A. 1.2 k 2 > 

If we allow the gluon to have a positive virtuality k 2 , (|A.3|) takes the form 



(*-*)"< (^-tt-jfl-j-^) (A.13) 
and leads to the following maximal invariant gluon mass for given Zi,z 2 : 

k 2 < k 2 m = (1 - zi)(l - z 2 ) - \ z x z 2 - 4m 2 - \j z x - Am 2 \J z 2 - 4m 2 

2 (1 - x + m 2 ) 



(?/+ — y)(y — y-) 

x — (2 — x)y — 4m 2 + */[ x 2 — 4m 2 ] [ (1 — y) 2 — Am 2 1 

(A.14) 

As we shall see below, the structure of the matrix element is such that essential contributions 
emerge from the logarithmic region y_ <3^y <3^y + , as well as from the vicinity of the lower limit, 
{y — y-) ~ y_~m 2 . Therefore for our purposes the following approximate expression suffices, 

k 2 m = (l-x)(y-y_) + 0(m 2 ) , (A.15) 



which differs from the exact formula ( [A. 14 ) only in a tiny region close to the upper integration 
limit, (y + — y) ~ m 2 . 

Two characteristic momentum scales emerge related to the limiting values of the ^-integration 

Q 2 (x) =W 2 -(l-x)y + « W 2 x(l-x); 

K 2 {x) =W 2 -{l-x)y^ « M 2 (l-x) 2 /x. 1 ' ' 

A.2 VASP Cross Sections 

The differential first order cross section integrated over angles of the final QQg system may be 
written as 

=<7-.%^IIc. (A.17) 
^ az\ az 2 J c 27r 

The subscript C marks the production current (vector (V), axial vector (A), scalar (S) or pseu- 
doscalar (P)) and the universal high-energy limit of the Born transition probability, 

Current (C) -> QQ. 

Straightforward calculation leads to the following expression for the II factors^: 



U c = 2( C S + H c ; (A. 18a) 

H V = (i=*+±Z*) = + (A .18b) 

\1— z 2 1— 2i/ y 1 - x 

Ha = (I + 2m 2 )H v + 4m 2 , H s = H P = H v + 2 . (A.18c) 



16 see pi for details 
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Here S is the "soft" bremsstrahlung term 
S 



_ (Vlii P2n\ 2 __ Zi + z 2 -l-2m 2 m 2 m 2 



Ui K 2 J (1- Zl )(l-Z 2 ) (1- Zl ) 2 (1-Z 2 ) 

x — 2m 2 1 m 2 1 — x + m 2 



x y y 2 (1 — x^ 2 



(A.19) 



Its contribution to each of the squared matrix elements is explicitly proportional to corresponding 
^-factor, the one that determines energy dependence of the corresponding QQ Born cross section, 



a, 



QQ 



^ ■ v ( c (v) ; 



Ca = Cs 



1 + 2m 2 , 
- 1 - 4m 2 



(A.20a) 
(A.20b) 







1 . 



As a result, the normalized differential distribution can be written in the following general form 
as a sum of a universal "iSoft" and a process dependent "7iard" contribution: 



d 2 wc = 



d 2 a 



c 



CfCIs 

2tt v 



{2S + (c 1 'Hc } dxdy. 



(A.21) 



z, G c Basis. Both structure of the matrix element and kinematics become particularly simple 
in terms of the gluon energy fraction and the gluon angle in the QQ cms (see ( |A.5| )). Making 
use of (|A.8|) one gets 



S 
7~Cv 



4(1 -z) (3 2 -u 2 



4(1-2) /3 2 sin 2 6 c 



4m 2 



4 



l-u 



2\2 



1 — U 2 



— 2 



4 



1 - (3 2 cos 2 e c 



(1 - /3 2 cos 2 G c ) 2 1 - ^ + 4m 2 ' 
2 . 



(A.22a) 
(A.22b) 



Invoking ( |A.7| ), ( A.21|) we obtain for the inclusive gluon energy spectrum in the vector channel 
(see also 12611) 



dw 



v 



C F a s dzp r 1 j 

/ ctcosBc < 2(1 — z) 



IT Z V J-l 



(3 2 sin 2 e c 



+ z' 



1 



1 



l-/? 2 cos 2 6 c 2 



Ck 1 . 



[l-(3 2 cos 2 9 c ) 2 

(A.23) 

Our convention to call the two pieces of the matrix element "soft" and "hard" becomes clear 
now: in the soft gluon limit, z<l, the second term is z 2 down compared to the first one. In 
this representation the "dead cone" phenomenon is also manifesting: the soft classical term S 
vanishes in the very forward directions (G c , it — @ c < 8q = m). 

In relativistic situation (m C ^) a collinear logarithmic enhancement occurs and the two 
pieces of ( |A.23j ) participate in forming the GLAP splitting function, 



del 

dw oc — - 



Co/-. \ 21 dz 1 + (1 

\2(l-z)+2r\ — oc — 

i j z 



dz . 
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It is worthwhile to notice that O (z" 1 ) and O (1) parts of the q — > q + g(z) splitting function 
are process independent, while the O (z) piece breaks factorization at the level of O (m 2 lnm 2 ) 
correction. Therefore the very notion of "fragmentation function" as a way of treating the jet 
evolution independently of the production mechanism gets lost beyond leading twist. 
To stress logarithmic character of the angular integration one may represent ( |A.23| ) as 



Cfol s dz r 1 dr) 



7TV z J no y/l - v 



2(l-z)^ + z 2 



1 1 

T] ~ 2 



c^ 1 



where 



r] = 1 - u 2 = 1 - (3 2 cos 2 6 C > 770 



4m 2 
1 - z 



(A.24a) 



(A.24b) 



For other production channels the second "hard" term in (|A.23|) , ( |A.24a| ) should be changed 
according to ( |Q8g ), ( |A.20b|) . 



A. 3 Integration over Virtual Boson Mass 

Considering inclusive characteristics {e.g., such as the quark energy spectrum) beyond the first 
order of perturbation theory one has to allow the gluon to decay in the final state, that is to 
have a positive virtuality k 2 , and to integrate over the latter within the available phase space. 
Actually, only such a combination of real and virtual gluon production leads to a sensible physical 
answer. First of all, an exclusive real gluon production cross section is clearly zero because of 
the standard infinite (double logarithmic) Sudakov form factor suppression. Secondly, and more 
importantly, a "real" gluon is an ill-defined object since its "on-mass-shell" interaction strength 
a may not be defined perturbatively (which is a real "infrared catastrophe" inherent for QCD). 

Integration over k 2 leads to appearance of the running coupling in the inclusive cross section. 
This nice property discovered in the pioneering papers by Gribov and Lipatov on partonic struc- 
ture of logarithmic field theories is particularly helpful in the QCD context. Here neither ao 
in "elastic" radiation nor small virtuality inelastic gluon decay systems are well defined. Mean- 
while, their combined action results in a legitimate a s (k m ) at the Euclidean scale related to the 
maximal available gluon virtuality. 

We shall demonstrate this correspondence taking care of subleading effects which might be 
relevant within the adopted approximation. To this end we write down a formal dispersion 
relation (with one subtraction) for the running coupling 

a(Q 2 ) = a ■ Z 3 (-Q 2 ) = a + Q 2 T ^ , (A.25a) 



h k 2 k 2 + Q 2 
where a is related to discontinuity of Z 3 at positive virtuality 



a{k 2 ) = ^ [Z 3 (k 2 - te) - Z 3 (k 2 + te)} = Disc {z 3 (k 2 )} . (A.25b) 



2ni 



CXq 



2tt 



In QED Z 3 is nothing but the photon renormalization function and cto ~ 1/137 — the on-mass- 
shell coupling constant. 
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Strictly speaking, such an identification is true for an Abelian theory only. In the QCD 
context the Ward identity between vertex and fermion propagator corrections, Z^Z q = 1, does 
not hold. As well known, both the non- Abelian vertex renormalization correction Z^ a ^ and the 
gluon propagator factor Z g participate (in a gauge dependent way) in forming the running a s . 
So one has to look upon Z% in the dispersion relation ( |A.25| ) as the gauge invariant product 



z 



(na) 



■ Zg ■ Z T 



(na) 



This subtlety does not affect the result, however. 

Another motivation is to use the nj dependence as a gauge to pinpoint the structure of the 
running a s in the inclusive cross section cJ. Quark loops (in order a 2 ) belong to Z g only, and 
employing a "quasi- Abelian" relation ( |A.25| ) becomes natural. 

In the problem under consideration the following structure emerges, 



/• fc m dk 2 o o 
a M( Zl , z 2 ; 0) + / — a(k 2 ) ■ M(z u z 2 ; k 2 ) 



a + 



k ldlf 



A- 2 



(A.26) 



a(k 2 ) }M( Zl ,z 2 ;0) + / dk 2 a(k 2 ) ■ AM(z u z 2 ; k' 



with k 2 m = k 2 n (z 1 ,z 2 ) the maximal squared virtual boson mass allowed for given x,y. Here we 
have singled out the k 2 = part of the matrix element by writing 

M( Zl , z 2 ; k 2 ) = M(z h z 2 - 0) + k 2 ■ AM(z u z 2 ; k 2 ) . (A.27) 

Now we may relate the characteristic integral in the first term in the rhs of (|A.26|) with the 
dispersion formula (|A.25a|) . To this end we employ the following exact representation of the 
dispersion relation 



a(Q 2 )=a + Q 2 



dk 2 a{k 2 ) 
o k 2 (k 2 + Q 2 ) 



a Q + 



Q 2 dk 2 , f2 , 



1 dt 



o 1 + t 



a(Q 2 /t)-a(Q 2 t) . (A.28) 



Perturbatively, a is of the order of a 2 . As a result, the last finite integral term constitutes a 
next-to-next-to-leading correction to the main one. To see this one has to view a as a slowly 
varying (logarithmic) function of its argument^ and to perform the Taylor expansion in In Q 2 to 
obtain 



i dt 



o 1 + t 



a{Q 2 /t) - a{Q 2 t) 



-2a'(Q 2 



o 1 + t o 



This exercise demonstrates a close correspondence between the spectral density ( |A.25b| ) and the 
/9-function @: 



a(Q 2 ) 

o-(Q 2 ) 



«o + J Q Q2 ^-°(k 2 ) + *-o\Q 2 ) + ... ; 
\Q 2 )-^ra"\Q 2 )+0{a^) . 



a 



(A.29a) 
(A.29b) 



17 which is true everywhere except in the very vicinity of a fermion threshold 
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Thus we express perturbatively the structure in curly brackets of (|A.26|) in terms of the running 
coupling at the Euclidean point Q 2 = k 2 



„2 . 



a(k 2 



a + 



a(k 2 



4m)-y«"(^) + 



(A.30) 



A. 4 Running Coupling in the First Order Spectrum 

It is straightforward to calculate the exact QQg matrix element with account of a non-zero virtual 
gluon mass k 2 > 0. ( |A.21| ) gets modified as follows: 



2S 
T~Cv 



- k 2 



HP + 2k 2 



1 



1 



(1-zi) 2 (l~z 2 ) 2 / 
zi + z 2 + k 2 



(l- Zl )(l- Z2 ) 



(A.31a) 
(A.31b) 



where and Tif are given by the original "on-mass-shell-boson" expressions ( |A.19| ) and 
( |A.18b|) respectively. In terms of ( |A.27p , 



A {25} 
&{H V } 



If 



(l-x) 2 1 



2(1 +x-y + k 2 
(l-x)y 



(A.32a) 
(A.32b) 



A. 4.1 k 2 = part of the matrix element. 



To evaluate the radiator (|2.7| ) we start by considering the first term in ( A.26Q proportional to 
M(z\, Z2', 0). One has to perform the y- integral with the factor a that emerges after integration 
over virtual boson mass, see (|A.30|) , 



y+ 



dya(k 2 m (x,y))-{2S^ + Cy 1 H^} 



(A.33) 



with k 2 m given by the approximate expression ( A. 15 ). Invoking ( A.19| ), ( |A.18b| ) we split the 
integrand into three pieces as follows: 



{2S^ + Cv 1 n {0) } 



x — 2m 2 
i 

1 — x 



+ G\l-x) 



2% 

y 

1 — x + m 2 x 
(1 — x) z 



1 — X 



(A.34a) 
(A.34b) 
(A.34c) 



41 



Logarithmic Piece ( |A.34a| ). This is the leading contribution to the radiator, in which the 
y- integration is logarithmic (collinear log): 



r ^ a(y - y_) . (A.35) 
Jy- y 



Here we have explicitly shown the essential y-dependence of the coupling factor a. The chain of 
approximations follows: 

(A.35) = r ^ L(y) + a'(y) In V -^- + ..Lf + * a(y) + a% C ^ ln(l - s) 

J y- y \ y j J y- y j t- s 



« r^a(2/)4«V) = / T a ^-ir a '^ + ^(«V) + « 3 (Q 2 )), 

jy_ ?/ 6 Jk 2 t 6 v 7 

(A.36) 

where ( |A.30|) has been used. We notice that a potential second order contribution a' oc a 2 at the 



low scale k 2 ~ M 2 cancels out, and the remaining term a'(Q 2 ) oc a 2 (W 2 ) is comparable with the 
two-loop correction to the hard cross section (coefficient function) and must be neglected within 
our approximation. 
Finally, 

f , . x fx — 2m 2 t . , \ rQ 2 dt , . 
y cfy (A34a) = f 2 — + C^(l - x) J y^ - a(t) . (A.37) 



Singular Piece ( [A.34bf ). The second term of ( |A.34j ) originates from the dead cone subtrac- 
tion. It is explicitly proportional to m 2 but this suppression gets compensated by the singular 
behaviour in y. Corresponding y-integral is concentrated in a tiny region (y — y-) ~ y- oc m 2 : 

r % My -v-) = «(».) % + Av-) r % "> — + ° <«"<»-)> < A - 38 > 

Jy- y Jy- y l Jy- V V- 

The O (a 2 (y-)) term vanishes in the relativistic approximation: 



^ln^= ^_*-ln a « f-A-ln S = 0. 
y- V (1 + s) 2 7o (1 + s) 2 

Thus, similarly to what has happened to the logarithmic term (|A.34aj) discussed above, a poten- 
tial O (a 2 (n)) contribution is absent. One is left with a pure a(n 2 ) correction to the radiator: 

dy (A346) = -2m 2 {—- — ) aU 2 ) = _ 2 ^ZM a r K *) . ( A .39) 

{ y- y+ J i - x 

Finite Piece ( [A.34cj ). Corresponding y-integral is collinear safe (i.e. finite in the m 2 = limit) 

and constitutes a(W 2 ) correction to the hard cross section. Here one extracts the coupling at 
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the upper limit, y ~ y + , and the ^/-integration becomes trivial (see ( |A.12| )): 



dy (A 34c) = a(Q 2 



~ z ~ — -+- C,y 



a 



[Q 2 )Vx 2 - Am 2 



[l-xf 
2 

1 — x 



+ c 



1 — x J 

2 ' 1-x 



_ x x — 2m 



2(1 — x) V 1 — x + 



v 



1 + 0(«(Q 2 

21 



(A.40) 



A. 4. 2 A part of the matrix element. 



From the first sight, corrections ( A. 32 ) to the QQg matrix element proportional to virtual boson 
mass look negligible: corresponding k 2 integration is no longer logarithmic, k 2 ~ k 2 m , and the 
result is of the order of cr(fc^) oc a 2 (/c^). 

This expectation indeed comes true for the "hard" correction term ( |A.32b| ) as well as for the 
second piece of (|A.32a|) . However, the first contribution to the "soft" correction term ( |A.32a|) 
is over-singular in y, as a result of which singularity logarithmic behaviour gets restored and a 
contribution to the two-loop anomalous dimension emerges. 

Indeed, 



fc m f /-fcm dk 2 k 2 

dkMk 2 ) = kiUki) + a'(ki)J o |rln-^ + 



k 2 m ■ a(k 2 J , (A.41a) 



"m "m 



where 

a = a' - a" + O (a"') . 
Invoking ( |A.15| ) for k 2 m one arrives at the following expression for the A-correction: 

fv+ dy 

- ( x - x ) / -*\v - y-) ■ °\v - v-) 

Jy- y 

An approximate evaluation follows: 

f y+d ^{y-y^)~a{y-y_)= f V+ %{y-y_)a{y) + 0{a\yJ)) 
J y- y J y- y 

V+ ^ a(y) - a{y_) = \a(Q 2 ) - a(n 2 ) - a'(Q 2 ) + a'( K 2 ) 
y- y 1 



(A.41b) 



(A.42) 



a'U 2 ) + 



(A.43) 



Again, as before, the a 2 (K 2 ) contribution cancels, and one finally gets the correction 
dy f dk 2 a(k 2 )AM(x,y;k 2 ) = -(1 - x) [a(Q 2 ) - a(K 2 ) ] + O (a 2 (Q 2 ) + cx 3 (k 2 ' 



(A.44) 



Combining (|A.37|) , (|A.39|) , ( |A.40| ) and (|A.44j) we finally arrive at the expression (|2.7|) for the 
perturbative radiator. 
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Appendix B : Second loop Anomalous Dimension A( 2 ) 



B.l AD (dimensional regularization) 

The two-loop anomalous dimension (AD) has been derived in the framework of dimensional 
regularization approach in fl9^pT||. In notations of Curci, Furmanski and Petronzio the non- 
singlet AD corresponding to Q — > Q transition in space-like evolution reads (eqs. (4.50) - (4.54) 

of H) 

1 + x 2 



where 



C F P qq (x, a) = aP(x) + a 2 7 



x 



P(x) 



X 



7 



W = C F P F (x) + C A P G (x)+n f T R P nf (x) 



(b.i; 



Pf 
Pg 
P n , 



-2P{x) In x ln(l — x) 



P(x) 

2 
~3 



1,2 11, 

i in x H 1k + 

z Q 1 



1 — x 

67 7T 2 

If 



2xj In x — |(1 + x) In 2 x — 5(1 — x) 
20 

f (1 + x) lnx H (1 — sc) 

3 



P(x) lnx+- +2(l-z 



(B.2a) 
(B.2b) 
(B.2c) 



Let us mention an extra contribution due to Q — > Q + QQ transition which reads 



Cpipgqix, a) = a 2 (C F - \C^j { 2P{-x)S 2 {x) + 2(1 + x) lnx + 4(1 - x) } , (B.3a) 



with 



So(x) 



(i+^O- 1 dz , 1 - 
— m 



(B.3b) 



This contribution (which formally belongs to the non-singlet AD, see, e.g., JT5]) vanishes as 
(1— x) 2 at large x, is colour suppressed and numerically negligible. In this paper it has been 
disregarded together with the singlet (sea) contribution to heavy quark yield. 

An algebraic massage leads to the following representation for the quark evolution kernef^j 



Cp Pqq(x, d) 



a + ICa 2 }P(x) + a 2 {a-C F V(x)+7Z(x)} - a' [P{x) \nx + 2(1 - ar) ] . (B.4) 



The first term here collects the one-loop AD with the part of 7^ 2 -* correction explicitly proportional 
to P(x) with the number /C given by ( |2.16b| ) above. These two combine forming a "physical" 
coupling in terms of the MS one according to ( |2.17b| ). This term totally absorbs the (1 — x)^ 1 
singularity of the evolution kernel: V(x) oc ln(l— x) and lZ(x) vanishes as (1—x) in the quasi-elastic 
limit. 



18 



The "+" prescription is implicit. 
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(The very last term in flB.4j) proportional to the derivative of the running coupling, 



a 



(I a s (Q) 



d\nW 2 2tt 



(tK ~ lT R n f ) a 2 + 



(B.5) 



does not count, as it is an artefact of the dimensional regularization scheme; see below.) 

The "true" second loop AD in curly brackets of ( |B.4| ) consists of two contributions. The first 
one is 



V{x) 



P(x) In 



x 



X 



^(1 + x) lnx 



lnx 



dz dy 5(x — yz) {P(y)}, P(z)\nz 
Jo Jo 



(B.6a) 



10 JO 

So defined, V has a very simple form in the moment representation, namely, 



V{j) = / dxx j ~ l V{x) 



P P 



Pi 



dxPix) 



dx 



3-1 



1 + x 2 

1 — X 



B.6b) 



(Obviously, [V(x)] + = V(x).) This (and only this) contribution changes (acquires an opposite 
sign) when time-like evolution is considered 0. Therefore ( |B.4| ) holds for both channels with 
a — ±1 referring to time- and space-like evolution correspondingly. The origin of the V term in 
the two-loop kernel may be traced backed to a simple kinematical difference between annihilation 
and scattering channels 111. It has been argued @ that a mismatch between the two-loop e + e~ 
and DIS anomalous dimensions would disappear (and thus the Gribov-Lipatov relation £1 would 
be rescued) if one considered scaling violation of "pseudo-moments" in x evaluated for a fixed 
value of {xW 2 } annih _ = {Q 2 Mscatt. 
Another term of the second loop AD is 

TZ(x) = (\C A -C F )P(x)\n 2 x-hC F [ \(l+x) lnx + (1-x) ] +C A [ (1+x) In a; + 3(l-x) ] . (B.7) 

Similarly to the one-loop splitting function P(x), it obeys the reciprocity relation 

- xTZix' 1 ) = K{x) , (B.8a) 
and, as a result, stays the same for time- and space-like evolution. Contrary to this, 

-xV(x~ 1 ) = [-l]-V(x). (B.8b) 



B.2 CF (dimensional regularization) 

O (a) correction to the hard cross section (jlB (coefficient function; CF) has to be taken into 
consideration together with the two-loop AD since neither CF nor AD is a scheme independent 
quantity beyond leading logs. 
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It reads (see eqs. (7.4), (7.10) of p[) 



+ 



8(1 — x) + aC F 
3(1 + xf 



T ± ^( ln V-^ + ^ (9 + 5x) 



1 — X 
aC F ■ 



\nx- |(1 + x) + vr 2 5(l - ar) 



(B.9a) 
(B.9b) 



The first line of ( |B.9aj ) is the CF for the scattering channel (Cf); the second line accounts for 
their difference. After simple manipulations one arrives at 



C A = (1 - %aC F )5(l-x) + aC F [P(x) (ln[x 2 (l - x)} - § ) - ±(5 + 9x) 



(B.10) 



In the cross section integrated over total QQg production angle (see ( |2.3j )) the combination 
C A +3C A emerges. With account of the longitudinal contribution the <5(1— x)-term gets modified, 



C A (x, a) = C A + 3Cf = (l-%aC F )5(l-x) + aC F [(B.10)]+ + 3aC F 



(1 + laC F )5{l-x) + aC F { « + 3 ; 



:b.id 



Constructing the CF that describes quark distribution normalized by the total cross section, 

C t {x, a) = (1 + ^clCf)- 1 ■ C A (x, a) , (B.12) 
one gets, with O (a) accuracy, the answer which may be presented in the following form: 



C t (x,a) = 6(l-x) + aC F [p(x) (ln[x(l - x)] - |) -±(l + x) 
+ aC F { P(x) lnx + 2(1 -x)} + . 

B.3 Scaling Violation Rate 

Scaling violation rate is an observable determined by the scheme invariant combination 



(B.13) 



D' 



d 



7 D(x,W 2 ) => P qq (x,a(W 2 )) + 



d 



dlnW 2 ' 

Combining (P-4j ) and ([B.13D one arrives at 



d\nW 2 



C t (x,a(W 2 )). 



(B.14) 



D' 



°MS + ^ a 



C F P{x) + a 2 CV { C F V{x) + ft(x) } 
+ a'C F [P(x) (ln[x(l-x)] 



(B.15) 



where the MS origin of the running coupling has been stressed in the relevant first order term. 



in 



The "+" prescription is implicit. 
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This result has to be compared with the scaling violation rate as described by ( p.7[ ). To this 
end one evaluates the In W 2 derivative of the relativistic radiator ( |2.10| ) to obtain 



D' =s- \a(Q 2 )C F P(x) - a'C F (l - x) + a 2 C F A^\ + a'C F {-^- + — - A . (B. 
L J [ 1 — x 2(1 — x) J 



16) 



In the leading term one has to expand the coupling of a composite argument Q 2 near W 2 as 
follows, 

a(Q 2 ) = a(W 2 ) + a' ln[x(l - x)] + . . . (B.17) 
Adding together terms proportional to a' one observes correspondence with a relevant part of 

(MM)- 

P(x) + (l-x), — r = iP(x)-i(l + x): 



l-x w v " 2(\-x) 4 w 4 

1 - x + 2(l-x 



In[x(l-x)] - (l-x) - + = P(x) (ln[x(l-x)] - f) - |(1 + s) 



Identifying 

(5.15) = (B.18) 

we arrive at the final result which relates an "effective" (dispersion scheme motivated) coupling 
a cB with a^jg and gives an expression for the "true" second loop correction to the radiator A^ 2 ) 
we were aiming at: 

a cS C F P(x) + a 2 C F A {2) = \a m + ICa 2 ] C F P(x) + a 2 C F { C F V(x) + K(x) } . (B.19) 
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